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MEASUREMENT  AND  INTERPRETATION  OF  NORTH  ATLANTIC 
OCEAN  MARINE  RADAR  SEA  SCATTER 


INTRODUCTION 

The  study  of  radar  sea  scatter  for  moderate  to  high  grazing  angles  has  received  much  attention 
over  the  past  several  years  because  of  airborne  and  spacebome  remote  sensing  and  military  tactical 
applications.  The  low  grazing  angle  case,  on  the  other  hand,  has  received  far  less  attention.  This 
case  is  more  complicated  for  several  reasons,  including  the  effects  of  shadowing  of  the  surface  by 
wave  crests;  the  occurrence  of  high-amplitude  transient  returns  referred  to  as  sea  spikes  [1-4];  and  the 
wider  range  of  sea  echo  strengths  expected  for  these  very  low  grazing  angles  as  a  result  of  modulation 
of  the  local  incidence  angle  by  the  slope  of  the  long  wave.  With  the  development  of  more  sophisti¬ 
cated  Navy  shipboard  radars,  it  has  become  apparent  that  a  satisfactory  sea  clutter  model  for  such 
radars  is  lacking,  thus  warranting  experimental  work  in  this  area.  In  addition  to  clutter  model  appli¬ 
cations,  the  shipboard  marine  navigation  radar  has  recently  received  attention  as  a  potential  tool  for 
remote  sensing  of  the  ocean  surface  by  using  its  surface-imaging  capabilities  [5,6]. 

In  the  work  reported  here,  statistics  of  the  normalized  radar  cross  section  (NRCS)  of  the  sea 
surface  were  determined  by  using  a  large  number  of  hourly  measurements  from  a  shipboard  marine 
navigation  radar.  The  data  were  collected  from  the  NOAA  research  vessel  Researcher  by  using  a 
horizontally  polarized,  high-resolution  X-band  radar  operating  at  grazing  angles  between  1°  and  7°. 
(For  this  report,  a  radar  is  considered  as  high  or  low  resolution  if  its  pulse  length  is  much  less  or 
much  greater  than  the  dominant  ocean  wavelength.)  The  amplitude  statistics  of  the  samples  of  the 
NRCS  of  the  ocean  surface  were  calculated  from  digitized  measurements  of  the  radar  video  signal  for 
each  of  a  series  of  range  bins  that  were  sorted  by  azimuth.  Data  from  a  60°  azimuthal  sector  were 
used,  the  sector  being  centered  about  the  direction  of  highest  radar  scatter  levels.  These  were  gen¬ 
erally  found  to  be  the  same  as  the  wind  direction.  Because  of  low  digital  recording  rates,  the  data 
were  collected  using  angle  increments  of  several  degrees  in  azimuth.  Many  antenna  rotations  were 
required  to  achieve  uniformly  filled  azimuthal  coverage.  Thus,  although  the  waves  are  not  imaged  as 
described  in  Ref.  5,  Statistics  of  the  scatter  from  all  portions  of  the  long  waves  were  sampled  as  the 
wave  system  traveled  through  a  range  bin  for  a  given  azimuthal  sector  over  many  antenna  rotations. 

Specifically  defined  parameters  of  the  NRCS  statistics  are  correlated  against  shipboard  measure¬ 
ments  of  wind  speed  and  air-sea  surface  temperature  difference.  From  these  quantities,  wind  stress  is 
calculated  and  compared  with  the  radar  measurements,  similar  to  the  procedure  used  for  airborne  and 
spacebome  scatterometers.  Applying  these  results  to  remote  sensing  of  wind  stress  leads  to  a  tech¬ 
nique  that  uses  a  relatively  low  data  rate  and  simple  signal  processing. 

EXPERIMENT  DESIGN 

Data  Acquisition 

An  experiment  was  designed  to  study  low  grazing  angle  sea  scatter  under  open  ocean  conditions 
under  joint  Navy  and  NOAA  cooperation.  The  NOAA  research  vessel  Researcher  was  instrumented  to 
collect  radar  sea  scatter  data  by  using  one  of  its  two  Raytheon  marine  navigation  radars,  the  antenna 
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of  which  is  shown  in  Fig.  1.  Although  the  Raytheon  radar  is  not  universally  used  on  oceangoing 
ships,  its  design  is  typical  of  such  radars,  and  it  is  worthwhile  to  discuss  its  operating  parameters  as 
they  affect  this  experiment.  The  shortest  pulse  available  from  this  system  is  65  ns  when  used  in  its 
short-range  mode.  This  mode  also  uses  a  3.5  kHz  pulse  repetition  frequency  (PRF),  which  yields  a 
286 -/is  period  between  pulses.  With  a  1°  antenna  beamwidth,  18  pulses  are  transmitted  while  the 
antenna  illuminates  a  target  with  a  1.8-s  rotation  period.  In  this  mode,  the  9.75-m  pulse  (65-ns  pulse 
width)  is  much  less  than  the  typical  open-ocean  dominant  wave  length,  and  the  radar  can  be  con¬ 
sidered  as  high  resolution  for  sea  scatter  studies,  according  to  the  definition  stated  previously. 


Fig.  1  —  The  radar  antenna  aboard  the  NOAA’s  ship  Researcher  is  a  typical 
marine  navigation  radar  antenna,  in  this  case  made  by  Raytheon 


A  data  acquisition  system  was  interfaced  to  the  output  of  the  Researcher’s  marine  navigation 
radars.  It  consisted  of  an  HP-86  computer,  a  Tektronix  7612D  transient  waveform  analyzer,  and  a 
Kennedy  tape  drive  with  a  Dylon  buffer  interface.  The  various  components  were  controlled  by  using 
an  IEEE-488  GP-IB  interface  system.  For  this  experiment,  data  processing  was  performed  offline  on 
a  VAX-750  computer. 

The  radar  was  operated  hourly  in  standard  fashion  on  a  not-to-interfere  basis,  except  for  selec¬ 
tion  of  the  minimum  range  display  (discussed  above)  while  recording  data.  The  data  were  digitized 
with  a  30-ns  spacing  over  256  samples  per  pulse,  for  a  high-resolution  measure  of  the  sea-surface 
echo.  These  were  recorded  with  an  azimuthal  spacing  of  7°,  determined  by  the  tape  drive  speed, 
with  3400  records  collected  over  60  rotations  of  the  antenna,  taking  about  2  min  of  each  hour. 
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The  system  was  calibrated  by  inserting  a  pulsed  X-band  signal  generator  output  into  the  direc¬ 
tional  coupler  at  the  radar  transceiver  box  and  incrementing  through  the  60  dB  dynamic  range  of  the 
receiver.  Absolute  calibration  was  then  achieved  by  towing  a  radar  calibration  sphere  (suspended 
from  a  weather  balloon  behind  a  boat)  through  the  antenna  beam  while  operating  the  antenna  in  a 
nonrotating  mode.  The  tow  line  used  was  sufficiently  long  to  allow  separation  of  the  boat’s  radar 
echo  from  the  calibration  sphere  as  the  two  passed  successively  through  the  antenna  beam. 

Data  Analysis  Technique 

The  first  stage  of  sea-scatter  data  processing  is  sorting  the  samples  according  to  antenna  pointing 
direction.  The  radar  sea  echo  maximizes  symmetrically  about  the  wind  direction,  and  this  echo  max¬ 
imum  was  used  to  select  the  region  of  scatter  to  be  analyzed.  Figure  2  shows  an  example  of  the 
range  of  received  power  encountered  vs  azimuth  for  one  depression  angle  or  range  bin  Data  from  a 
60°  sector  centered  about  the  azimuth  of  wind-generated  maximum  scatter  were  chosen  for  statistical 
analysis  for  each  of  10  range  bin  samples.  These  bins  geometrically  correspond  to  10  depression 
angles,  because  the  antenna  pattern  had  essentially  constant  gain  with  elevation  angle  over  the  angles 
illuminated.  This  design  maintains  a  uniform  radar  PPI  display  while  the  ship  is  pitching  and  rolling 
in  high  seas.  Received-power  data  from  this  chosen  azimuthal  sector  were  first  converted  to  normal¬ 
ized  radar  cross  section  (NRCS)  for  each  data  point  by  using  receiver  calibrations  and  the  measured 
cross  section  of  the  calibrated  sphere.  Next,  a  cumulative  histogram  of  NRCS  was  assembled  tor 
each  range  bin  and  plotted  on  Rayleigh  probability  paper.  On  such  plotting  paper,  the  Rayleigh  dis¬ 
tribution  plots  as  a  straight  line  with  unit  slope.  The  Weibull  distributions  plot  linearly  also,  but  the 
slope  is  a  variable  parameter  defined  by  the  distribution.  The  Weibull  distribution  yields  a  realistic  fit 
for  sea  scatter  [7-9]  and  allows  parameterization  of  the  NRCS  distribution  by  using  two  variables. 
The  cumulative  probability  density  function  of  linear  NRCS  values  x  is  represented  by  the  Weibull 
distribution,  parametric  in  A  and  B,  and  is  given  by 

p(x)  =  exp  {-(xfl)  /  A}.  (1) 

(See,  for  example.  Ref.  10,  p.  109.)  The  data  are  plotted  as  NRCS  in  dB  vs  the  expected  value, 
10*  log  (jc).  From  Eq.  (1),  this  is  log  [  -  In  [pOO]  -  log  A .  The  Rayleigh  probability  axis  is  gen¬ 
erated  by  substituting  fractional  values  between,  but  not  equal  to,  0  and  1  (e.g.,  0.999  and  0.001). 
On  such  Rayleigh  paper,  the  slope  of  the  straight  line  is  B,  which  is  also  called  the  shape  parameter. 
The  parameter  A  is  called  the  scale  parameter.  Both  can  be  related  to  the  mean  value  of  the  distribu¬ 
tion  n  by  fi  =  A~bT(B  +  1),  where  T  is  the  T  function. 

The  distribution  function  of  NRCS  values 

dp  (x)  dx  =  (B /A)xb~x  exp  [  ~(xB)  I A  )  (2) 


is  the  derivative  of  Eq.  (1). 

Figure  3(a)  plots  two  theoretical  Weibull  distributions  using  Eq.  (2)  and  the  parameters  listed, 
each  with  a  different  percentage  scaling  factor.  Figure  3(b)  shows  a  composite  plot  of  these  two  dis¬ 
tributions,  consisting  of  one  of  the  distributions  in  one  percentage  region  and  the  second  distribution 
in  the  second  region.  The  cutoff  between  the  two  distributions  is  obvious.  The  distributions  were 
generated  on  a  linear  NRCS  basis,  but  they  are  plotted  vs  the  log  of  NRCS,  as  with  real  data.  This 
plot  is  shown  to  illustrate  the  reason  for  plotting  the  cumulative  distribution  rather  than  the  probability 
distribution  for  the  data.  When  plotting  the  probability  distribution,  it  is  not  a  simple  matter  to  dif¬ 
ferentiate  between  the  two  distributions  for  two  reasons.  First,  the  higher  amplitude  contribution 
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Fig.  2  —  The  amplitude  distribution  of  received  power  of  low  grazing 
angle  sea  scatter  is  shown  as  a  function  of  azimuthal  angle.  This  peak  in 
azimuth  is  probably  determined  by  the  direction  of  the  large  scale  waves, 
typically  coincident  with  the  wind  direction,  but  it  can  lag  the  wind 
direction  under  changing  conditions. 
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Fig.  3(a)  —  A  pair  of  Weibull  distributions  generated 
from  Eq.  (4)  by  using  the  Weibull  parameters 
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Fig.  3(b)  —  A  single  distribution  created  from  using 
each  of  the  distributions  of  Fie.  3(a).  With  such  a 
distribution  of  real  data,  the  break  in  the  curve 
would  be  difficult  to  identify  using  the  PDF  due  to 
inherent  noisiness. 
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generally  consists  of  a  few  to  40%  of  the  data,  and  the  peak  location  is  simply  not  present  in  the  com¬ 
posite  plot.  Second,  because  of  the  general  noisiness  of  real  data,  the  break  point  of  transition  of  one 
curve  to  the  other  (evident  in  Fig.  3(b))  is  not  obvious  either.  However,  because  of  the  markedly  dif¬ 
ferent  tails  of  the  two  distributions,  the  cumulative  plot  on  Rayleigh  paper  allows  this  differentiation 
between  Weibull  distributions  to  be  made  relatively  easily,  with  as  little  as  5%  of  the  data  points  in 
the  highest  amplitude  contribution.  Note  that  simply  adding  the  two  distributions  does  not  provide  a 
pair  of  straight  lines  on  the  Rayleigh  paper;  it  provides  a  curve  with  a  hyperbolic  shape  near  the  tran¬ 
sition  from  one  Weibull  to  the  other.  To  achieve  two  intersecting  straight  lines,  as  evidenced  in  the 
data,  the  first  distribution  must  be  held  up  to  the  level  of  intersection  of  the  two  curves,  at  which 
point  the  second  distribution  is  in  effect.  It  is  as  though  one  type  of  scattering  feature  simply  replaces 
the  first  over  an  appropriate  percentage  of  the  surface,  affecting  the  NRCS  distributions  accordingly. 

Figure  4  shows  the  result  of  plotting  three  typical  cumulative  distributions  of  NRCS  values  for  a 
single  elevation  angle  (radar  time  delay)  for  low  to  moderate  winds  on  Rayleigh  paper.  For  the 
second  and  third  case,  there  are  three  different  regions  of  the  plotted  data,  the  highest  two  of  which 
are  fitted  by  straight  lines.  For  very  low  wind  conditions,  the  regions  to  either  side  of  the  middle 
linear  region  disappear,  resulting  in  a  single  Weibull  distribution. 
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Fig.  4  —  A  cumulative  probability  distribution  of  sea  surface  NRCS 
data  points.  The  resulting  function  appears  bimodal.  with  the  high 
NRCS  region  identified  as  sea  spikes  and  the  middle  NRCS  region 
identified  with  a  distributed  scatter  source,  such  as  Bragg  scatter. 
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Software  was  developed  for  least-squares  fitting  of  straight  lines  to  both  linear  regions  of  the 
cumulative  distribution,  beginning  with  the  high  NRCS  region.  The  linear  fits  to  these  two  regions 
were  defined  for  postanalysis  by  six  parameters:  (a)  a  slope  (or  shape  parameter)  B  and  (b)  a  50% 
intercept  for  the  lower  amplitude  line;  (c)  the  NRCS  amplitude  and  (d)  percentage  value  of  the  inter¬ 
section  of  the  two  fitted  straight  lines;  and  (e)  the  slope  B,  and  (f)  the  “effective  mean  value”  of 
NRCS  of  the  higher  amplitude  straight-line  fit.  This  effective  mean  value  is  taken  as  the  decibel- 
mean  of  the  NRCSs  at  the  two-line  intersection  and  the  1  %  intersection  of  the  high-amplitude-region 
linear  fit.  These  results  provide  a  complete  statistical  characterization  of  the  data  and  are  stored  for 
later  correlation  with  surface  truth  parameters. 

General  Behavior  of  the  Data 

Before  comparing  radar  and  environmental  data,  some  of  the  characteristics  of  the  data  should 
be  interpreted  in  light  of  possible  scattering  mechanisms.  The  higher  amplitude  contribution  in  Fig.  4 
are  associated  with  so-called  radar  sea  spikes,  which  have  been  identified  in  previous  experiments 
[1-4].  In  these  experiments,  the  radar  was  not  rotating  but  illuminated  the  same  patch  of  sea  surface 
continuously.  The  term  sea  spike  is  used  because  of  the  temporal  behavior  of  certain  periods  of  the 
radar  return.  During  these  periods,  the  signal  is  much  stronger  and  has  a  short-term  amplitude  varia¬ 
tion.  This  variation  is  much  less  than  the  scatter  signal  that  varies  in  a  periodic  manner  and  is  similar 
to  that  of  the  sea  surface.  Thus  sea  spikes  are  expected  to  occupy  the  highest  amplitude  region  of  the 
probability  distribution. 

Although  they  have  been  studied  for  shoaling  conditions  near  shore  by  Kalmykov  and  Pustovoy- 
tenko  [11]  and  the  NRL  group  [1-4],  except  for  the  work  of  Hansen  and  Cavaleri  in  the  Gulf  of 
Mexico  [12],  few  deep-water  experiments  have  been  conducted.  The  observations  of  Hansen  and 
Cavaleri  have  deterministically  associated  high-amplitude  radar  returns  with  white-water  features 
around  tne  crests  uf  waves.  These  observations  have  been  made  by  using  split-screen  visual  presenta¬ 
tions  of  the  radar  cell  and  the  amplitude  time  sequence  of  the  received  echo  power.  In  Weibull  fitting 
of  their  data,  they  also  showed  a  break  in  the  high-amplitude  region  but  did  not  deterministically  asso¬ 
ciate  individual  contributions  to  the  distribution  with  specific  echos  in  time.  We  assume  that  the 
higher  amplitude  distribution  in  our  data  is  also  caused  by  crest  features.  This  assumption  is  verified 
by  the  behavior  of  the  data  with  increasing  scattering  area. 

The  behavior  of  the  data  in  the  lowest  amplitude  region  of  Fig.  4  looks  similar  to  the  tail  of  a 
Gaussian  distribution  when  plotted  on  Rayleigh  paper.  We  suggest  that  this  region  contains  samples 
of  receiver  noise  that  result  for  data  samples  recorded  at  time  delays  occurring  in  the  shadowed  area 
of  the  radar  illumination,  just  behind  the  wave  crest.  This  is  discussed  in  Appendix  A.  This  feature 
is  found  to  disappear  for  very  low  winds  and  waves,  as  would  be  expected  as  shadowing  diminishes. 

The  middle  region  for  depression  angles  >2°  generally  encompasses  the  majority  of  the  data 
samples.  This  region  may  be  associated  with  the  major  fraction  of  the  illuminated  surface,  where  a 
composite-scatter  model  would  be  expected  to  hold  when  a  distributed  scattering  mechanism  is  dom¬ 
inant.  This  region  reduces  to  the  single  distribution  found  for  mild  winds  when  shadowing  and  wave 
breaking  are  minimal.  Further  evidence  for  these  preliminary  interpretations  follows. 

RESULTS 

For  correlation  with  the  radar  data,  environmental  measurements  were  made  of  temperatures  of 
the  air  and  seawater,  as  well  as  wind  speed  and  direction.  Wind  speed  was  measured  by  an  anemom¬ 
eter  at  a  height  of  30  m  on  the  ship’s  forward  mast.  Wet  and  dry  bulb  air  temperatures  were  mea¬ 
sured  hourly  on  the  bridge.  Seawater  temperature  was  measured  by  a  water  intake  port  at  the  bow  of 
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the  ship  at  a  depth  of  -2  m.  Data  were  collected  throughout  a  6-day  period  in  a  region  ranging  from 
long.  40°  to  42°  W  and  lat.  71°  to  72°  N,  off  the  coast  of  Long  Island,  New  York.  Results  are 
shown  in  Fig.  5. 
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Fig.  5  —  A  plot  of  the  environmental  parameters  measured  aboard  the  NOAA  ship 
Researcher  for  the  days  covered  by  the  experiment.  The  days  covering  the  radar 
data  show  a  relatively  constant  wind  direction  throughout  the  period. 


During  the  6  days,  the  first  3-day  period  featured  winds  blowing  from  2  to  18  m/s  from  the 
northwest,  off  the  coast  of  Long  Island,  resulting  in  a  minimum  fetch  of  185  km.  On  the  fourth  day, 
the  winds  changed  direction  by  180°  as  a  result  of  the  passage  of  a  front.  This  gave  radar  results  that 
differed  markedly  from  the  first  3-day  period  because  of  confused  seas.  On  10  February,  additional 
low-wind-speed  data  were  available  for  analysis  as  well,  before  a  slow  and  continuous  change 
occurred  in  wind  direction.  The  radar  data  for  the  period  of  changing  winds  show  a  great  deal  of 
scatter,  representing  a  nonequilibrium  condition  between  winds  and  waves.  We  consider  here  only 
the  cases  for  near-equilibrium  conditions,  using  data  collected  during  the  period  between  7  February 
and  1900  local  time  on  9  February,  and  for  a  short  period  on  10  February. 

The  reduced  NRCS  distribution  parameters  described  earlier  are  first  compared  with  wind 
speed.  The  data  are  differentiated  according  to  air-sea  temperature  differences,  with  different  plotting 
symbols  used  for  2°  increments  in  temperature  [13].  Results  for  the  various  NRCS  parameters  varia¬ 
tion  with  wind  speed  follow. 

NRCS  Median  vs  Wind  Speed 

Consider  first  the  variation  of  the  50%  intercept,  or  median,  of  the  cumulative  distribution  for 
individual  incidence  angles.  (Note  that  this  median  generally  is  also  the  median  of  the  midlevel  distri¬ 
bution,  unless  the  high-level  region  exceeds  50%. )  The  median  is  plotted  vs  wind  speed  in  Fig.  6  for 
four  of  the  10  processed  range  bins  corresponding  to  depression  angles  of  1.0°,  1.3°,  1.6°,  2.0°, 
2.4°,  3.1°,  3.8°,  5.1°,  6.6°,  and  7.7°.  (These  angles  were  chosen  for  analysis  to  provide  a  loga¬ 
rithmic  variation  in  depression  angle  spacing.)  Each  set  of  points  has  been  least-squares  fit  to  a 
straight  line.  The  slope,  variance,  and  10-m/s  wind  speed  intercept  is  shown  at  the  bottom  of  each 
figure.  (Note  that  this  linear  fitting  corresponds  to  conventional  treatment  of  depression-angle  varia 
tion  of  NRCS,  both  in  the  area  of  scatterometry  and  low-grazing-angle  sea-scatter  modeling.  It 
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presupposes  a  power-law  dependence  of  NRCS  on  wind  speed,  an  assumption  being  reconsidered  in 
some  circles).  A  comparison  with  previous  measurements  of  Sittrop  |14|.  who  used  a  pulse  length  of 
120  ns  and  HH  polarization  shows  excellent  agreement  with  our  results.  This  pulse  length  is  roughly 
twice  the  value  of  ours  but  is  still  in  the  high-spatial-resolution  regime,  as  defined  earlier. 

The  last  two  depression  angles,  6.5°  and  7.7°,  correspond  to  time  delays  that  are  affected  by  the 
time-dependent  receiver  attenuation  known  as  sensitivity-time-control  (STC),  which  is  used  to  protect 
the  receiver  from  leak-through  of  the  transmitted  pulse.  Although  the  slope  of  median  NRCS  vs  wind 
speed  is  valid  for  these  two  depression  angles,  the  absolute  values  should  be  scaled  upward  by  an 
unknown  amount,  since  receiver  calibrations  are  performed  at  longer  time  delays. 

Depression  Angle  Dependence  of  10-m/s  Fitted  Value 

The  results  of  the  linear  fits  of  Fig.  6  have  been  summarized  in  Fig.  7  in  the  form  of  the  depen¬ 
dence  of  the  intercepts  and  slopes  on  the  depression  angle.  For  the  10-m/s  intercepts  of  NRCS,  note 
the  very  weak  depression  angle  dependence,  with  just  four  decibels  separating  the  first  value  from  the 
eighth.  This  agrees  with  the  distributed  scatter  arguments  made  earlier.  The  last  two  angles  now 
show  more  quantitatively  the  effects  of  the  STC  attenuation  discussed  earlier.  The  error  bars  on  each 
point  are  the  standard  errors  of  the  linear  fits.  The  results  of  Fig.  7(a)  are  fitted  by  eye  to  a  straight 
line,  with  a  0.45  slope  on  the  log-log  plot. 

Figure  7(a)  shows  that  the  NRCS  values  occupy  similar  amplitude  levels  from  one  grazing  angle 
to  the  next  and  exhibit  similar  slopes  at  the  lowest  grazing  angles.  This  apparent  similar  NRCS  varia¬ 
tion  over  these  five  lowest  depression  angles  can  be  described  by  one  of  two  different  models.  Either 
the  median  NRCS  is  independent  of  both  illuminated  area  and  grazing  angle  or  the  NRCS  variation 
with  each  of  the  two  parameters  causes  their  mutual  cancel!  ttion,  probably  a  less  likely  phenomenon. 
The  Bragg-scatter  basis  of  the  composite  scatter  model  predicts  area  independence  of  the  NRCS. 
However,  the  NRCS  values  and  the  weak  incidence-angle  dependence  is  definitely  not  in  agreement 
with  the  evidence  discussed  in  Appendix  A  for  horizontal  polarization.  (Again  with  the  caveat  that 
medians  for  high-spatial-resolution  data  and  means  of  low-spatial-resolution  data  are  not  necessarily 
comparable.)  The  very  weak  variation  of  NRCS  with  increasing  radar  cell  size  gives  evidence  for  a 
distributed  surface  scattering  mechanism  in  the  midamplitude  region.  Therefore  we  refer  to  this 
scatter  as  the  distributed  scatter  region. 

Figure  7(b)  shows  a  plot  of  the  slopes  of  Fig.  6  vs  wind  speed.  These  results  for  the  last  two 
angles  are  valid  now.  since  the  change  of  the  median  NRCS  with  wind  speed  is  independent  of  STC 
attenuation.  It  also  shows  that  the  shadowing  is  not  affecting  the  median  NRCS-wind  speed  slope  vs 
depression  angle. 

Empirical  Model 

For  engineering  design  interests,  an  empirical  model  can  be  formulated  to  summarize  the  results 
of  the  last  two  figures.  A  straight  line  has  been  fitted  to  the  data  of  Fig.  7: 

log  [NRCS|Om/,,(0)!  =  -3.18  +  0.45  log  [0(deg)],  (3a) 


m(0)  =  2.5  +  5.0  log  (0(deg)) 


1  <  0  <  4.5 


(3b) 
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Pig,  6  —  The  50%  intercept,  or  median  value,  of  the  hourly  radar  Matter  measurement!  plotted  vi  wind  speed 
for  4  of  10  radar  depression  angles  for  the  period  covering  2/6  (2 100)  to  2/10  (2100) 
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Fig.  7(a)  —  Distributed  scatter  NRCS  value  at  the  10- 
m/s  intercept  plotted  vs  radar  depression  angle.  The 
last  two  points  show  the  effects  of  the  radar  STC 
circuit  at  the  shortest  ranges  sampled.  The  two  curves 
were  calculated  using  the  composite  scatter  model  with 
tilt  modulation,  adding  a  mean  surface  slope  of  16° 
and  20°  to  each  depression  angle  labeled  on  the  x  -axis. 


Fig.  7(b)  —  The  slope  of  the  linear  fits  of  Fig.  8 
plotted  vs  radar  depression  angle,  indicating  the 
wind  speed  variation  of  the  NRCS  change  for  each 
angle 


The  variation  of  the  median  NRCS  of  the  distributed  scatter  is  thus  given  by;. 

log  (NRCS^/NRCS i on,/,)  =  m  (0)*  log  (Uw/ 10)  (4) 

where  Uw  is  the  wind  speed  measured  aboard  ship. 

Percentage  Value  of  Sea  Spikes 

The  next  parameter  to  consider  is  what  may  be  viewed  as  the  percentage  occurrence  of  sea 
spikes,  defined  by  the  value  of  the  intersection  of  the  two  fitted  lines  on  the  Rayleigh  paper,  indicated 
by  the  dotted  lines  in  Fig.  4.  This  percentage  value  can  also  be  considered  as  the  fractional  percen¬ 
tage  of  NRCS  values  associated  with  the  higher  amplitude  linear-fit  region  of  the  data.  Results  for 
the  10  depression  angles  presented  above  are  plotted  vs  the  log  of  wind  speed  (Fig.  8).  This  time,  a 
line  of  3.75  slope  has  been  drawn  through  the  data,  based  on  a  model  for  wave-breaking  occurrence 
to  be  discussed  later.  No  attempt  has  been  made  to  optimize  the  fits  by  choosing  an  appropriate  ori¬ 
gin  for  each  case,  and  the  intercept  at  100%  occurs  at  ~  10- m/s  wind  speed. 

It  is  apparent  that  the  data  for  the  higher  grazing  angles  more  closely  fit  a  3.75  slope  than  those 
for  lower  grazing  angles.  This  is  expected  because  the  percentage  occurrence  of  sea  spikes  is  a  func¬ 
tion  of  both  the  occurrence  of  scatterers  on  the  surface  and  the  occurrence  of  shadowing.  Shadowing 
is  expected  to  have  a  stronger  effect  at  the  lowest  depression  angles  for  high  waves.  That  is,  for  the 
higher  wind  speeds,  sharper  crested  waves  will  produce  more  shadowing  at  the  lowest  grazing  angles, 
thus  reducing  the  observed  percentage  occurrence  of  sea  spikes. 
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Fig.  8  —  Plots  of  the  percentage  coordinate  of  the  intersection  of  the  two  linear  fits,  interpreted  as  the  percentage  occurrence 
of  sea  spikes,  plotted  vs  wind  speed  for  the  10  radar  depression  angles.  For  angles  <2°,  the  high-wind-speed  results  fall 
below  the  dotted  line,  perhaps  as  the  result  of  shadowing,  as  discussed  in  the  text. 
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Fig.  8  —  (Continued)  Plots  of  the  percentage  coordinate  of  th.'  intersection  of  the  two  linear  fits,  interpreted  as  the 
percentage  occurrence  of  sea  spikes,  plotted  vs  wind  speed  for  thj  10  radar  depression  angles.  For  angles  <2°,  the  high- 
wind-speed  results  fall  below  the  dotted  line,  perhaps  as  the  result  of  shadowing,  as  discussed  in  the  text. 
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Fig.  8  —  (Continued)  Plots  of  the  percentage  coordinate  of  the  intersection  of  the  two  linear  fits,  interpreted  as  the 
percentage  occurrence  of  sea  spikes,  plotted  vs  wind  speed  for  the  10  radar  depression  angles.  For  angles  <2°,  the  high- 
wihd-speed  results  fall  below  the  dotted  line,  perhaps  as  the  result  of  shadowing,  as  discussed  in  the  text. 


Although  there  is  much  scatter  in  the  data  and  not  all  air-sea  temperature  ranges  are  available 
for  all  wind  speeds,  a  weak  trend  in  temperature  effects  is  discernible,  but  the  expected  scatter  in 
wind  speed  precludes  drawing  any  firm  conclusions  for  this  parameter.  The  wind  speed  error  is  dis¬ 
cussed  in  greater  detail  later. 

A  Model  for  U*75  Wind  Speed  Dependence 

Consider  the  apparent  C/*‘75  dependence  of  the  percent-intercept  of  Fig.  8.  Hansen  and  Cavaleri 
[12]  provide  evidence  that  radar  sea  spikes  are  associated  with  features  near  the  crests  of  waves,  par¬ 
ticularly  those  crests  that  are  nearly  or  actually  breaking.  Whitecap  coverage  is  also  associated  with 
breaking  waves  and  the  residual  white  water.  Thus,  the  percentage  occurrence  of  whitecap  coverage 
and  the  percentage  occurrence  of  radar  sea  spikes  are  expected  to  scale  in  a  similar  manner  because 
they  are  both  related  to  breaking  or  nearly  breaking  waves.  As  discussed  in  Appendix  B,  Wu  [18] 
derived  a  model  predicting  a  3.75-power  wind-speed  dependence  for  whitecap  coverage  data,  which  is 
the  reason  for  the  3.75-slope  line  plotted  over  the  sea  spike  data  reported  here.  Thus,  the  reason  for 
this  3.75-slope  fit  is  a  theoretical  one,  implicitly  associating  percentage  of  radar  sea  spikes  with  per¬ 
centage  whitecap  coverage  as  events  with  a  common  physical  origin  but  without  attempting  to  relate 
the  proportionality  factor  between  the  two  at  this  time.  Further  radar  experiments  are  planned  to 
determine  whitecap  coverage. 
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Effects  of  Shadowing  of  Small  Crests 

For  grazing  angles  >2°,  the  fit  of  the  data  to  the  line  is  good.  For  angles  of  <2°,  either 
shadowing  of  the  surface  is  affecting  the  fit  for  the  higher  wind  speeds  or  a  different  scattering  or 
propagation  phenomenon  is  affecting  the  data.  Shadowing  could  be  expected  to  cause  the  observed 
reduction  of  sea  spike  percentages  at  high  winds  and  steep  waves  for  the  lowest  depression  angles.  It 
would  do  this  by  reducing  the  radar  illumination  of  smaller  crests  on  the  steep  forward  face  of  the 
dominant  wave  that  occur  as  a  result  of  breaking  (the  spilling  breakers  considered  by  Wetzel  [20]). 
(Appendix  A  describes  this  low-grazing-angle  sea-scattering  background.  Figure  A  1(c)  shows  the  sea 
surface  illuminated  for  a  very  low  grazing  angle.)  For  this  extremely  low  depression  angle,  only  the 
crests  of  the  dominant  waves  are  being  illuminated.  For  weaker  winds,  less  developed  seas,  and  the 
same  depression  angle  shown  in  Fig.  Al(c),  the  highest  crests  would  not  be  as  steep.  This  would 
cause  less  shadowing,  and  a  higher  fraction  of  crests  present  on  the  forward  face  of  the  longest  waves 
would  be  illuminated.  Thus,  for  a  given  radar  depression  angle,  the  fraction  of  wave  crests  on  the 
surface  that  are  not  shadowed  is  dependent  on  long  wave  steepness,  and  implicitly,  on  the  wind  speed 
and  time  which  it  has  blown. 

For  Figs.  8(e)  through  8(j),  apparently  all  crests  are  still  illuminated  for  the  wind  speeds 
reported  and  shadowing  reduces  illumination  only  of  the  distributed  scattering  surface  area.  This 
increases  the  percentage  value  at  which  the  top  of  the  cumulative  distribution  of  Fig.  3  rounds  off  to 
noise  samples  but  only  at  the  expense  of  number  of  samples  of  distributed  scatter.  At  2°,  it  appears 
as  though  there  is  a  transition  region  to  shadowing  of  smaller  wave  crests,  since  the  data  for  the  lower 
wind  speeds  lie  reasonably  close  to  the  3.75-slope  line,  and  the  waves  for  these  winds  still  allow  most 
of  the  crests  present  to  be  illuminated. 

Predictions  of  scattering  effects  occurring  near  2°  have  been  published  in  the  literature.  Radar 
propagation  effects  near  this  angle  have  also  been  observed,  such  as  the  observed  change  in  range 
falloff  of  NRCS  from  R-4  to  R  8  (see,  for  example,  Long  [17]).  This  effect  has  previously  been 
suggested  as  being  the  result  of  interference  effects  between  the  forward  scattered  wave  and  the  direct 
wave,  assuming  geometric  optics  holds.  Wetzel  has  proposed  a  shadowing  function  for  low-grazing- 
angle  scatter  to  explain  such  propagation  effects  [16].  The  model  suggests  effective  “islands  of 
scatter”  as  a  result  of  large  crest  shadowing.  These  islands  are  just  the  peaks  of  the  dominant  waves 
in  his  simplified  formalism.  He  has  suggested  that  this  mechanism  is  responsible  for  the  change  in 
range  dependence  near  the  2°  depression  angle.  The  only  assumption  Wetzel  makes  about  the  sea 
surface  is  one  of  mean  slope  appropriate  to  a  fully  developed  sea.  We  suggest  here  that  the  percen¬ 
tage  sea  spike  falling  below  the  dashed  lines  for  the  lowest  grazing  angles  in  Fig.  8  is  also  a  result  of 
the  onset  of  shadowing  of  smaller  wave  crests  by  those  associated  with  the  dominant  waves,  another 
manifestation  of  the  shadowing  model  proposed  by  Wetzel. 

Reconciliation  of  U0  5  and  Ul  Dependences  of  Drag  Coefficients 

Beyond  this  implicit  relationship  of  radar  sea  spikes  to  whitecap  coverage,  one  can  go  further  in 
suggesting  a  reconciliation  of  the  i/J  dependence,  typically  observed  for  drag  coefficient  variation,  vs 
the  Uw5  required  by  Wu.  (Uw5  is  necessary  for  the  U^  75  dependence  of  whitecap  coverage  developed 
in  Appendix  B,  which  also  fits  our  sea  spike  results  reasonably  well.)  In  Appendix  B,  we  review 
Donelan’s  suggestions  [19]  that  the  total  rate  of  energy  exchange  between  winds  and  waves  can  be 
broken  up  into  a  sum  of  several  contributions,  one  of  which  is  wave  breaking.  Whereas  the  total  drag 
coefficient  measured  varies  linearly  with  U ,  we  suggest  that  the  wave  breaking  portion  varies  as  U  \ 
and  this  contribution  is  used  for  the  radar  sea  spike  and  whitecap  measurements.  Thus,  although 
Wu’s  derivation  is  based  on  the  total  wind  stress,  only  the  contribution  due  to  wave  breaking  should 
be  used  in  the  derivation.  In  fact,  it  may  be  that  the  large  spread  of  data  points  typically  found  in 
determining  a  drag  coefficient,  by  empirically  plotting  the  wind  friction  velocity  vs  wind  speed,  is  due 
to  the  variation  in  the  contributions  of  Eq.  (B5)  for  the  experimental  conditions  encountered. 
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Sea  Spike  RCS  vs  Wind  Speed 

Effective  Mean  of  the  Sea  Spike  Region  vs  Wind  Speed 

The  dB-mean  of  the  sea  spike  1%  intercept  and  the  intercept  of  the  two  linear  fits  can  be  also 
plotted  vs  wind  speed.  Figure  9  shows  these  results  for  four  of  the  incidence  angles.  For  angles 
>2°,  the  data  vary  quite  strongly  with  incidence  angle,  in  contrast  to  the  distributed  scatter  results. 

Mean  Sea  Spike  NRCS  at  10  m/s 

This  incidence  angle  variation  is  the  10-m/s  intercept  values,  shown  plotted  vs  depression  angle 
in  Fig.  10(a).  The  sharp  change  is  seen  at  2°.  Dotted  circles  have  been  drawn  for  the  last  two 
angles.  These  angles  use  corrections  of  4  and  8  dB,  the  same  values  used  earlier  for  the  distributed 
scatter  case.  Figure  10(b)  shows  these  same  data  converted  to  absolute  cross  section.  At  <2°,  the 
RCS  is  nearly  constant  near  10  m2,  rising  sharply  to  a  peak  near  2.5°  and  dropping  precipitously  with 
increasing  angle.  Finally,  Fig.  10(c)  shows  the  change  in  sea  spike  NRCS  with  wind  speed,  also 
determined  from  Fig.  9  and  the  other  available  angles. 

1  %  Sea  Spike  Intercept 

To  give  an  indication,  for  radar  design  interests,  of  the  largest  sea  spike  to  be  expected  as  a 
function  of  wind  speed,  the  1%  intercepts  of  the  distributions  are  shown  plotted  vs  wind  speed  as 
absolute  RCS  in  dBm2  in  Fig.  11.  The  results  for  all  10  grazing  angles  are  further  summarized  in 
Fig.  12  as  10-m/s  values  of  RCS  plotted  vs  depression  angle.  At  <2°,  cross  sections  of  up  to  30  m2 
(15  dBm2)  can  be  expected  1%  of  the  time  for  10-m/s  winds.  At  >2°,  the  same  discontinuity 
observed  before,  the  maximum  values  grow  quite  markedly.  Absolute  cross  sections  of  the  order  of 
1000  m2  are  seen  to  occur  for  6°  and  7°  depression  angles.  This  is  not  important  for  shipboard 
operation  because  the  ranges  are  just  a  few  hundred  meters,  but  it  could  be  important  for  airborne  and 
space-based  radars. 

Sea  Spike  Model  for  RCS  Depression  Angle  Characteristics 

A  model  is  proposed  to  explain  the  above  results  that  is  a  variation  on  one  first  published  by 
Wetzel  [20].  The  wedge  model  of  Lyzenga  et  al.  [21]  yields  RCS  that  increases  with  increasing 
depression  angle  for  horizontal  (HH)  polarization,  rather  than  the  falloff  above  4°  that  we  observe, 
and  it  does  not  describe  the  sharp  change  near  2°.  Wetzel’s  model  has  difficulty  with  this  discon¬ 
tinuity  as  well,  but  his  model  can  be  modified  to  account  for  it. 

Figure  13  shows  the  geometry  of  three  models  for  the  plume  of  water  spilling  down  the  front 
face  of  a  long  wave,  the  residual  body  of  water  after  a  crest  has  broken.  0  is  the  angle  that  the  toe  of 
the  plume  makes  with  the  surface,  a  is  the  local  slope  of  the  long  wave,  and  0g  is  the  local  grazing 
angle.  The  model  shown  in  Fig.  13(a)  is  that  chosen  by  Wetzel  and  has  a  value  of  0  equal  to  zero. 
This  choice  for  0  forces  the  toe  to  be  perpendicular  to  the  local  long  wave  at  the  intersection.  This 
assumption  allows  scatter  to  the  plume  from  the  ground  plane  on  the  near  flat  surface  forward  of  the 
plane.  This  creates  a  comer  reflector  geometry,  which  in  turn  allows  a  variation  of  RCS  with  depres¬ 
sion  angle  to  occur.  Scatter  can  occur  for  all  incidence  angles,  from  a  surface  wave  parallel  to  the 
ground  plane  with  negative  incidence  angle  on  upward.  However,  such  a  plume  model  is  not  stable 
because  of  downward  gravitational  forces,  and  one  would  expect  a  collapse.  This  collapse  would 
create  a  toe  that  would  be,  at  a  minimum,  perpendicular  to  the  vertical,  as  shown  in  the  second 
model,  Fig.  13(b).  For  this  case,  0  —  a,  and  a  comer  reflector  mode  is  not  allowed,  i.e.,  no  scatter 
can  occur  from  the  fiat  surface  ahead  of  the  plume,  then  off  the  plume,  then  back  to  the  radar.  (The 
comer  reflector  mode  disappears  for  all  0  >  0).  For  the  choice  of  0  made  here,  with  the  toe  perpen¬ 
dicular  to  the  horizontal,  scatter  can  occur  only  for  positive  grazing  angles.  Neither  of  these  two 
models  predicts  the  RCS  discontinuity  that  we  observe  near  2°. 
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Fig.  9  —  Plots  of  the  "effective  mean"  of  the  sea  spike  distribution,  as  defined  in  the  text,  plotted  vs  wind  speed 
for  four  radar  depression  angles.  Note  the  definite  dependence  on  depression  angle. 
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Fig.  10(a)  —  The  sea  spike  mean  10-m/s  intercept 
of  Fig.  9  plotted  vs  radar  depression  angle.  The 
last  two  points  show  the  effects  of  the  radar  STC 
circuit,  with  dotted  data  points  also  drawn  in, 
applying  the  correction  estimated  from  Fig.  7(a). 
A  sharp  change  is  seen  in  the  region  around  2°, 
similar  to  previous  parameter  changes  in  this 
region. 


Fig.  10(b)  —  The  sea  spike  total  RCS  value  at  the 
lOm/s  intercept,  results  of  Fig.  10(a)  multiplied 
times  the  scattering  cell  size,  plotted  vs  radar 
depression  angle.  The  rightmost  two  points  have 
been  corrected  and  fit  quite  well  the  negative  2.33 
slope  with  depression  angle.  The  sharp  change  is 
still  seen  in  the  region  around  2°. 


2  3  4  56789  10 

DEPRESSION  ANGLE  (deg) 


Fig.  10(c)  —  The  change  in  sea  spike  NRCS 
with  wind  speed  plotted  vs  depression  angle 
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Fig.  11  —  The  lee  spike  absolute  RCS  at  the  1%  intercept  of  the  Weibull  distribution  plotted  vs  wind  speed  for  four  of  the 
10  depression  angles  used,  giving  an  estimate  of  the  clutter  level  attained  1%  of  the  time  for  detection  statistics 
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Fig.  13(b)  —  Modified  plume  model,  with  angle  of  plume  intersection 
with  ground  plane  perpendicular  to  the  horizontal.  This  model  does  not 
allow  scatter  from  the  ground  plane  forward  of  the  plume  but  only  from 
the  specular  line  along  the  cylinder  normal  to  the  incoming  ray. 
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To  predict  such  a  discontinuity,  one  must  assume  that  the  toe  of  the  plume  is  such  that  its  nor¬ 
mal  is  2°  positive  to  the  horizontal,  regardless  of  the  slope  of  the  long  wave  as  shown  in  Fig.  13(c). 
Energy  is  reflected  directly  back  to  the  radar  only  for  specular  reflection  from  the  cylinder  face,  i.e.. 
only  for  rays  whose  extensions  pass  through  the  center  of  curvature  of  the  cylinder.  Hence,  for  the  (3 
chosen,  this  occurs  only  for  angles  greater  than  the  0^,0  shown  in  the  figure.  The  toe  angle  (3  must 
therefore  be  equal  to  the  sum  of  the  slope  of  the  long  wave  at  the  toe  a  plus  the  minimum  grazing 
angle  which  is  2°  in  our  case. 

This  model  now  predicts  a  constant  cross  section  with  increasing  grazing  angle,  since  reflections 
from  the  long  wave  ground  plane  do  not  provide  a  scattering  mode.  The  expression  for  such  scatter 
from  a  cylinder  is  given  by 


RCSCV|  =  kaL 2 


(4) 


where  a  is  the  radius  of  curvature  of  the  cylinder,  k  is  the  radar  wave  number,  and  L  is  the  median 
length  of  a  plume.  If  we  set  the  RCS  equal  to  63  m2  (18  dBm2),  the  following  combinations  of  a  and 
L  satisfy  the  equation  at  X-band  for  mean  RCS  values  above: 

a  (cm)  13.2[6.4|  12.8 
L( m)  1 3.2 1 2.2 1  1.6 

Using  the  1000-m2  result  for  the  1%  intercept,  the  following  crest  lengths  result: 

a  (cm)  |  3.2|6.4|  12.8 
L( m)  |  12.6 1  8.9 1  6.3' 

These  are  well  within  realistic  surface  features  for  a  breaking  plume,  and  all  values  of  a  satisfy 
assumptions  of  geometrical  optics. 

Although  we  have  produced  realistic  plume  features  to  describe  the  maximum  RCS  observed, 
the  RCS  is  constant  for  grazing  angles  above  9g  with  this  model.  However,  we  can  account  for  the 
sharp  falloff  beyond  2°  if  we  consider  the  radar  range  at  which  this  occurs,  the  beamwidth  of  the 
radar,  and  the  scale  size  of  the  plume  length  across  the  radar  beam.  For  a  10-m  crest  length  and  a  1  ° 
azimuthal  beamwidth,  the  crest  fills  the  beam  at  a  range  of  573  m.  For  the  antenna  height  of  24  m, 
this  corresponds  to  a  grazing  angle  of  2.4°.  For  shorter  ranges,  the  crest  more  than  fills  the  beam 
and  is  no  longer  fully  illuminated,  now  acting  as  a  distributed  scatter  in  the  cross  range  dimension. 
The  illuminated  portion  is  now  the  cross  range  dimension: 
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In  addition,  we  can  write 


R  =  H  / sin  4> 


(6) 


in  the  flat  earth  approximation  (<5%  error  for  angles  >2.5°),  and 

R  ~  H /</> 


(7) 
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in  the  small  angle  approximation.  The  expected  decrease  in  cross  section  with  decreasing  range  as 
the  result  of  incomplete  crest  illumination  is  therefore  given  by: 

a  =  (2ira/\)(R  d<t>)2.  (8) 

When  range  is  converted  to  depression  angle,  this  becomes 

a  =  (2xa/\)H2  d<t>2  (9) 

and  partially  accounts  for  the  observed  angular  falloff  above  2.5°  (-2  of  the  -2.33  power  observed). 
By  using  either  Wetzel’s  model  or  that  of  Lyzenga  et  al.,  increasing  RCS  is  predicted  through  this 
region,  which  is  in  disagreement  with  our  observations. 

Finally,  for  grazing  angles  <2°,  our  model  predicts  no  scatter  from  the  plume.  However,  sea 
spikes  are  still  present  according  to  the  observed  double  Weibull  distribution,  albeit  with  the  discon¬ 
tinuous  drop  in  RCS.  Recall  that  these  Weibull  distributions  are  a  result  of  sampling  over  the  surface 
profile  with  a  spatial  scale  smaller  than  the  dominant  wave.  We  have  identified  the  first  Weibull  dis¬ 
tribution  with  distributed  scatter,  presumably  from  portions  of  the  wave  away  from  the  crest.  The 
RCS  distribution  is  determined  primarily  by  the  long  wave  slope  distribution  according  to  our  compo¬ 
site  model.  In  similar  fashion,  we  can  associate  the  second  distribution  with  crest  regions:  with 
plumes  providing  the  RCS  at  angles  >2°  and  a  weaker  crest  scattering  mechanism  for  angles  <2° 
(actually  for  all  angles,  but  dominated  by  the  plume  model  at  angles  >2°). 

For  this  region,  we  consider  the  wedge  model  proposed  by  Lyzenga  et  al.,  who  calculated  the 
cross  section  per  unit  crest  length-squared,  RCS/L2.  This  expression  evaluated  at  0°  depression  angle 
is  -29,  -23.  and  -17  dB  for  external  wedge  angles  of  195°,  210°,  and  240°,  respectively.  For  the 
10-m  crest  length  used,  these  become  -9,  -3,  and  +3  dBm2,  respectively.  Our  mean  values  lie  at 
~  12.5  dBm2,  while  the  1%  intercept  values  are  at  15  dBm2.  Thus  a  wave  crest  length  longer  than 
10  m  is  necessary  to  satisfactorily  describe  our  results  when  using  wedge  scatter.  However,  the 
numbers  used  earlier  as  lengths  are  plume  lengths,  which  can  be  expected  to  be  less  than  the  crest 
length  that  created  them,  according  to  Wetzel’s  development.  With  this  caveat,  we  have  produced  a 
model  that  satisfactorily  interprets  all  the  characteristics  of  our  observed  sea  spike  scatter  amplitudes. 

Distributed-Scatter  Weibull  Shape  Parameters 

A  second  parameter  required  to  fully  describe  the  Weibull  parameterization  of  the  data,  in  addi¬ 
tion  to  the  median  values,  is  the  shape  parameter  B  of  Eq.  (1).  This  shape  parameter  is  the  slope  of 
the  linear  fit  of  Fig.  4.  It  can  be  expected  to  vary  from  unity,  appropriate  to  a  normal  distribution  of 
radar  echo  field  strengths  (expected  to  occur  for  very  light  winds),  to  our  maximum  observed  values 
of  the  order  of  4.  This  parameter  that  increases  beyond  cr.c  with  larger  wind  speed  indicates  the 
increasing  dynamic  range  of  NRCS  over  which  the  distribution  varies.  Although  the  parameter 
against  which  the  slope  is  plotted  is  the  wind  speed,  long-wave  slope  increasing  with  wind  speed  is 
probably  the  cause  of  the  growth  in  slope.  This  is  expected  according  to  the  two-scale  model  because 
the  local  incidence  angle  at  the  scattering  patch  is  increased  on  the  average,  which  in  turn  contributes 
a  larger  NRCS.  (This  is  more  fully  discussed  in  Appendix  A.) 

Figure  14  shows  results  for  this  parameter  for  one  of  the  10  depression  angles.  For  very  low 
grazing  angles,  the  shape  parameter  is  of  the  order  of  unity  with  little  wind  speed  dependence.  While 
the  median  NRCS  increases  with  wind  speed,  the  NRCS  spread  or  dynamic  range  does  not  do  so  for 
low  depression  angles.  It  is  possible  that  this  is  an  effect  of  shadowing  of  the  major  portion  of  the 
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Fig.  14  -  The  distributed  scatter  Weibull  shape  parameter  B ,  determined  from  the  slope  of  the  linear  fit  of  Fig.  5, 

plotted  vs  wind  speed  for  four  of  the  10  depression  angles 
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sea  surface  by  the  wave  crests.  Less  of  the  wave  surface  is  illuminated,  and  as  a  result,  the  selection 
of  local  incidence  angles  is  limited.  These  results  can  be  summarized  again  by  plotting  their  10-m's 
intercepts  and  slopes  vs  depression  angle.  No  models  currently  exist  to  explain  this  behavior,  since 
no  derivation  of  the  Weibull  distribution  or  any  other  distribution  from  first  principles  has  been 
accomplished.  These  results  provide  a  first  comparison  set  for  such  a  model. 

Figure  15(a)  plots  the  10  m/s  intercepts  vs  depression  angle.  Again,  a  sharp  change  is  seen  near 
2°.  For  all  angles  >2°,  the  value  of  B  at  10  m/s  wind  speed  is  very  nearly  3.  Perhaps  this  is  fortui¬ 
tous,  since  the  slope,  or  rate  of  increase  of  B  with  wind  speed,  varies  with  grazing  angle  (Fig.  15(b)). 
At  <2°,  we  have  fit  both  a  linear  fit  as  well  as  one  that  asymptotically  approaches  unity  for  small 
grazing  angles,  i.e.,  a  Rayleigh  distribution. 
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Fig.  15(a)  —  Intercepts  at  10  m/s  of  the  distributed  Fig.  15(b)  —  The  change  in  Weibull  slope  parameter 

scatter  Weibull  shape  parameter  plotted  vs  depressions  with  wind  speed  plotted  vs  depression  angle 

angle,  indicating  the  wind  speed  variation  of  B  for 
each  angle 

Sea  Spike  Scatter  Weibull  Shape  Parameters 

Figure  16  shows  an  example  of  the  data  for  sea  spike  Weibull  slopes.  The  10-m/s  intercept  of 
these  fits  to  the  Weibull  B  parameter  is  shown  in  Fig.  17(a)  vs  depression  angle.  Again,  a  sharp 
change  is  seen  near  2°.  For  higher  angles,  the  shape  parameter  is  remarkably  constant  about  the 
value  2.  The  slope  of  the  Weibull  parameter— wind  speed  variation  is  shown  in  Fig.  17(b)  to  com¬ 
plete  the  summary.  Again,  no  models  are  available  that  explain  this  behavior. 

Median  NRCS  vs  Friction  Velocity  for  Four  Drag  Coefficients 

To  account  for  the  differences  in  atmospheric  stability  that  occurred  during  the  experiment,  the 
measured  air-sea  temperature  differences  were  used  in  four  different  models  to  provide  a 
temperature-dependent  drag  coefficient.  This  attempt  to  account  for  atmospheric  stability  in  air-sea 
interaction  studies  is  driven  by  remote  sensing  needs,  and  the  models  that  we  use  have  all  been 
developed  in  conjunction  with  radar  experiments.  These  areas  of  research  include  scatterometry  and 
synthetic  aperture  radar  image  interpretation. 
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DEPRESSION  ANGLE  (deg)  DEPRESSION  ANGLE  (deg) 

Fig.  17(a)  —  Intercepts  at  10  m/s  of  the  sea  spike  Fig.  17(b)  —  Slope  of  the  sea  spike  Weibull  shape 

Weibull  shape  parameter  plotted  vs  depression  angle.  parameter  plotted  vs  depression  angle,  indicating  the 

At  >2°  the  parameter  is  nearly  constant,  as  was  wind  speed  variation  of  B  for  each  angle.  These 

observed  with  the  distributed  scatter  Weibull  parameters  are  very  small  at  <2°. 

parameter.  At  <2°,  the  shape  parameters  approach 
one,  indicating  a  Rayleigh  distribution. 

The  drag  coefficient  allows  the  radar  data  to  be  parameterized  against  wind  friction  velocity. 
This  is  done  by  the  bulk  method  of  estimating  the  friction  velocity  that  uses  empirically  measured 
drag  coefficients.  The  bulk  method  assumes  an  equilibrium  between  winds  and  waves  in  terms  of 
momentum  transfer  between  them.  For  rapidly  changing  conditions  such  as  at  the  onset  or  reversal  of 
winds,  more  accurate  techniques  are  required.  The  dissipation  technique  measures  the  fluctuations  of 
the  wind  field,  while  the  eddy  correlation  technique  measures  three  orthogonal  components  of  the 
wind.  Because  each  of  these  techniques  requires  a  very  stable  platform,  they  are  typically  not  suitable 
for  shipboard  operation. 

The  first  temperature  drag  coefficient  used  was  that  of  Weissman  et  al.  [22];  the  second,  that  of 
Keller,  Plant,  and  Weissman  (KPW)  [13]  but  using  the  V 0  5  wind  speed-dependent  drag  coefficient  of 
Wu  discussed  in  Appendix  B;  the  third,  that  of  KPW  using  the  drag  coefficient  of  Large  and  Pond 
[37];  and  the  fourth,  that  of  Wu  [23].  Figure  18  shows  results  for  a  depression  angie  of  1.02.  Fig. 
19  shows  results  for  a  depression  angle  of  7.7°  as  plots  of  median  NRCS  vs  wind  friction  velocity. 
The  wind  friction  velocities  were  determined  by  using  the  measured  wind  speeds  and  Eq.  (B3)  of 
Appendix  B. 

The  greater  air-sea  temperature  difference  increases  the  wind  friction  velocity  compared  to  that 
value  for  zero  air-sea  temperature  difference.  Thus,  for  the  data  analyzed,  the  effects  are  only  to 
spread  the  data  in  the  horizontal  plane  because  there  is  a  strong  correlation  between  wind  speed  and 
air-sea  temperature  difference.  It  is  not  clear  whether  the  apparent  )0-cm/s  limit  in  the  data  for  7.68° 
is  real,  since  there  is  a  large  spread  in  the  data  points  in  dB  for  the  2°  to  0°  range.  The  data  at  the 
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Fig.  18  —  Median  NRCS  value  at  1.02°  depression  angle  vs  wind  friction  velocity,  for  four  different 
drag-coefficient  models  that  account  for  atmospheric  stability  effects 
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Fig.  19  —  Median  NRCS  value  at  7.68°  depression  angle  vs  wind  friction  velocity,  for  four  different 
drag-efficient  models  that  account  for  atmospheric  stability  effects 
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lowest  wind  speeds  were  collected  on  10  February,  after  the  change  in  wind  direction.  Perhaps  this 
large  spread  is  the  results  of  ambient  waves  from  the  previous  wind  field  propagating  into  the  area. 

Generally,  conclusions  cannot  be  drawn  as  to  whether  there  is  an  air-sea  temperature  difference 
effect  from  these  data  alone.  Nevertheless,  they  provide  preliminary  empirical  data  on  which  to  base 
low  grazing  angle  scatterometry  when  using  a  marine  radar. 

RELATED  EXPERIMENTS 

This  section  discusses  several  radar  sea  scatter  papers  that  contain  pertinent  results.  When  con¬ 
sidered  as  a  group,  results  of  these  experiments  allow  inferences  that  can  be  used  to  form  the  basis  of 
a  general  radar  scattering  model  and  that  explain  some  of  the  features  of  our  data.  These  experiments 
include  short-pulse  airborne  Doppler  spectrum  measurements  [24];  dual-polarized,  short-pulse  non¬ 
coherent  measurements  [12];  modeling  of  long-pulse,  coherent  radar  Doppler  spectral  features  as  a 
function  of  sea-surface  parameters  [25];  dual-polarized  Doppler  spectra  for  wind-driven  waves  in  a 
wave  tan’  [26];  and  comparison  of  horizontal  HH  and  vertical  VV  polarized  Doppler  spectra  near  a 
coastline  for  relatively  low  wind  conditions  [27]. 

Airborne  Coherent  Doppler  Measurements 

Hicks  et  al.  [24]  conducted  an  airborne  radar  experiment  using  a  high  incidence  angle,  high  spa¬ 
tial  resolution,  and  horizontal  polarization.  They  found  that  returns  from  regions  near  the  crests  of 
waves  gave  a  larger  radar  return  and  a  larger  Doppler  shift  than  returns  from  points  along  the  rest  of 
the  wave.  Although  these  results  are  for  higher  grazing  angles  than  those  considered  here,  they  are 
the  only  deep-water  results  that  both  identify  the  strongest  radar  return  with  breaking  features  on  the 
waves  and  associate  a  maximum  Doppler  shift  with  that  region  of  the  wave.  Because  the  sampling 
scheme  of  our  data  does  not  identify  the  source  of  the  large-NRCS  region  of  scatter,  we  can  rely  on 
the  observations  of  Hicks  to  identify  this  region  with  crests. 

Dual-Polarized  Measurements 

The  dual-polarized  noncoherent  radar  sea  scatter  measurements  of  Hansen  and  Cavaleri  [12] 
were  nonsimultaneous  but  were  conducted  near  enough  in  time  so  that  the  environmental  conditions 
did  not  change  significantly  from  one  polarization  to  the  next.  The  data  were  analyzed  in  a  manner 
similar  to  that  used  here.  Cumulative  distributions  were  calculated  for  the  two  polarizations,  and  an 
example  of  their  results  is  shown  in  Fig.  20.  Note  that  the  sense  of  the  slope  change  is  the  result  of 
the  method  of  parameterization— plotting  the  percentage  above  a  given  value  rather  than  below  the 
value.  The  breaks  in  our  data  are  found  to  occur  at  the  same  percentage  and  NRCS  point  when  plot¬ 
ted  either  way.  However,  it  was  found  that  the  break  between  the  two  distributions  could  be  identi¬ 
fied  by  curve- fitting  software  more  effectively  when  plotted  by  using  the  technique  in  our  report,  par¬ 
ticularly  for  low  percentage  occurrence  of  sea  spikes.  This  is  because  the  percentage  of  sea  spikes 
are  spread  over  a  wider  range  in  the  horizontal  dimension  when  plotted  our  way. 

The  important  point  to  notice  in  this  figure  is  that  throughout  the  region  of  moderate  NRCS 
values  or  Bragg  region  the  distributions  for  the  two  polarizations  have  similar  slopes  but  different 
50%  intercepts  or  medians.  This  consistent  difference  in  NRCS  throughout  this  region  are  what 
might  be  expected  from  Bragg  scatter  predictions  for  HH  and  VV  polarizations  in  the  low-grazing- 
angle  regime.  In  addition,  if  the  slope  of  the  Weibull  fit  is  determined  primarily  from  tilt  modulation 
by  long  waves,  the  statistics  of  all  W  echoes  and  the  HH  distributed  scatter  echoes  should  be  deter¬ 
mined  by  the  statistics  of  the  long  wave  slope  distribution,  resulting  in  similar  Weibull  slopes  for  both 
polarizations.  This  picture  is  very  much  in  agreement  with  these  observations. 
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(a) 


Fig.  20  —  Cumulative  distributions  of  NRCS  are  shown  for  HH  and  VV 
polarizations  for  X-band  radar  scatter  collected  by  Hansen  and  Cavaleri  in  the 
Gulf  of  Mexico.  Note  the  similar  slopes  in  the  moderate  NRCS  region  of 
scatter,  in  agreement  with  Eq.  (3). 


The  second  point  to  note  about  the  figure  is  that,  while  the  HH  data  show  a  break  to  a  second 
Weibull  fit  in  the  high-amplitude  sea-spike  region  as  our  results  do,  the  VV  results  show  no  such 
break.  The  VV  data  in  the  high  NRCS  region  of  the  distribution  continues  to  fit  the  line  defined  by 
the  moderate  NRCS  values.  Since  the  same  ensemble  of  surface  features  are  illuminated  for  either 
polarization,  differences  in  the  highest  amplitude  NRCS  values  for  the  two  polarizations  indicate  a 
different  scattering  mechanism  for  VV  polarization  than  for  HH  polarization  by  the  surface  features 
that  are  responsible  for  the  HH  sea  spikes.  Hansen  and  Caveleri  found  that  the  sea  spike  returns  as 
seen  on  a  time  history  of  radar  echo  corresponded  to  times  when  wave  crests  passed  through  the 
scattering  cell,  in  agreement  with  observations  of  Hicks  et  al.  (Thus  the  crests  apparently  scatter 
HH-polarized  energy  differently  than  VV  polarization  and  behave  as  distributed  scatterers  for  the  VV 
polarization,  but  more  coherent  ones  for  HH  polarization.) 
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Doppler  Spectrum  Peak  Shift  Model 

Reference  25  reanalyzes  X-band  data  taken  from  a  published  Russian  experiment  [28].  In  Ref. 
28,  Doppler  shifts  for  HH,  HV,  and  VV  polarizations  are  measured,  together  with  surface  truth  con¬ 
sisting  of  wind  speed,  rms  wave  height,  and  dominant  wave  period.  The  reanalysis  found  that  for 
VV  polarization  the  Doppler  shift  of  the  peak  of  the  Doppler  spectrum  reported  could  be  explained  by 
a  two-scale  model.  In  this  two-scale  model,  the  scatterers’  Doppler  shifts  were  determined  primarily 
by  the  orbital  motions  of  the  large  waves  and  by  the  various  surface  currents  that  could  arise,  includ¬ 
ing  Stokes  currents  caused  by  the  nonsinusoidal  profile  that  occurs  for  higher  waves.  These  results 
are  in  contrast  to  the  conclusions  drawn  by  the  Russian  authors  that  the  composite  scatter  model  could 
not  be  used  to  explain  their  measurement  results.  However,  they  had  not  considered  a  sufficient 
number  of  effects  that  could  contribute  to  the  Doppler  shift  imposed  on  the  radar  echo  by  the  sea  sur¬ 
face.  For  light  winds,  where  its  contribution  is  a  significant  fraction  of  the  total,  it  is  necessary  to 
include  the  phase  velocity  of  the  Bragg- resonant  capillary  waves  to  fit  the  data.  This  provides 
evidence  for  a  model  of  Bragg  scatter  from  wind-driven  gravity  capillary  waves. 

For  HH  polarization,  higher  scatterer  velocities  observed  cannot  be  explained  by  orbital  waves 
and  currents,  but  the  superimposed  motions  of  accelerating  wave  crests  on  the  large  waves  should  be 
considered.  The  HH  Doppler  shifts  scale  linearly  when  plotted  vs  the  VV  shifts,  with  a  constant  mul¬ 
tiplier  of  1.4.  Crest  features  responsible  for  the  HH  scatter  are  moving  faster  than  the  orbital  wave 
velocity,  as  would  be  expected  for  a  breaking  wave.  This  agrees  with  the  results  of  Hicks  et  al.,  in 
their  observation  of  the  highest  Doppler  shifts  and  largest  echo  power  levels  from  the  region  of  the 
wave  crest  for  HH  polarization.  The  types  of  crest  motions  required  to  satisfy  the  1 .4  scaling  factor 
observed  between  HH  and  VV  Doppler  shifts  were  not  considered. 

Doppler  Spectra 

The  results  discussed  so  far  are  similar  to  the  double-peaked  Doppler  spectra  of  Kwoh  and  Lake 
[26]  for  HH  and  VV  data.  These  data  cover  a  time  containing  the  passage  of  several  long  paddle- 
driven  waves  with  wind,  produced  in  a  wave  tank.  Figure  21  shows  Doppler  spectra  for  VV  and  HH 
polarizations,  respectively.  Both  spectra  exhibit  double  peaks,  but  the  high  and  low  peaks  are 
exchanged  for  horizontal  and  vertical  polarization.  The  magnitude  of  the  shift  of  the  HH  peak  is 
larger  than  the  1 .4  scaling  factor  presented  by  Melnichuk  and  Chemikov.  However,  the  wave  spectra 
produced  in  a  wave  tank  are  known  to  be  quite  different  from  open-sea  spectra,  so  such  a  comparison 
is  probably  not  quantitatively  meaningful.  Nevertheless,  the  sense  of  the  magnitude  of  the  shift  of  the 
dominant  peak  in  the  spectrum  is  the  same  as  observed  in  Ref.  28,  indicating  similar  scattering 
mechanisms  in  each  case.  Since  the  spectra  of  Kwoh  and  Lake  included  the  passage  of  several  long 
waves,  a  Doppler  shift  could  not  be  assigned  to  a  specific  portion  of  the  wave  as  was  done  in  Hicks, 
et  al. 

Open  Ocean  X-Band  Doppler  Spectra 

Figure  22  shows  spectra  collected  by  Shibata  et  al.  [27]  with  a  C-band  coherent  radar  and  a 
0  5-fis  pulse  for  very  mild  winds  of  4  m/s.  Approach-and-recede  contributions  to  the  spectrum  are 
observed  that  are  unambiguous  in  Doppler  frequency.  These  occur  in  a  manner  similar  to  those  typi¬ 
cally  observed  for  HF  radar  spectra.  The  peaks  at  ±150  Hz  appear  as  minor  contributions  to  the 
Doppler  spectrum  but  are  at  -1.4  times  the  peak  of  the  main  contribution.  This  agrees  with  the 
HH/VV  Doppler  shift  ratio  of  Ref.  28,  but  the  larger  Doppler-shifted  peak  is  not  the  dominant  peak 
in  the  spectrum,  as  in  Ref.  28.  However,  all  results  discussed  so  far  have  been  for  X-band  radars,  in 
which  case  sharp  crest  features  may  contribute  a  greater  amplitude  than  for  lower  radar  frequencies. 
This  spectrum  is  also  observed  to  be  quite  different  from  the  wave  tank  measurements  of  Kwoh  and 
Lake,  which  were  also  made  at  X-band. 
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Fig.  21  —  Radar  Doppler  spectra  for  scatter  from 
waves  in  a  wave  tank  for  HH  and  VV  polarizations. 
Note  the  bimodal  character  of  each  spectrum  with 
similar  location  of  peaks  but  of  differing  amplitudes. 
(Figure  3  from  ‘‘The  Nature  of  Microwave  Back- 
scattering  From  Water  Waves,"  by  D.S.  Kwoh.  et 
al.,  from  The  Ocean  Surface,  p.  254,  published  by 
D.  Reidel  Publishing  Company,  1985,  used  by 
premission.  | 
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Fig.  22  —  Radar  Doppler  spectra  for  scatter  from  open  ocean  waves  for  a  C- 
band  radar.  The  peaks  in  Doppler  scale  in  a  similar  fashion  as  reported  by 
Melnichuk  and  Chemikov  at  X-band  for  HH  and  VV  polarization.  [Figure  1 
from  “Doppler  Spectra  of  Microwave  Radar  Echo  Returned  From  Calm  and 
Rough  Sea  Surfaces,”  by  A.  Shibata,  et.  al.,  from  The  Ocean  Surface,  p.  265, 
published  by  D.  Reidel  Publishing  Company,  Copywright  1985,  used  by  prem¬ 
ission.] 


Summary  of  Experiments— Properties  of  a  Doppler  Shift  Model 

Considered  as  a  group,  these  results  show  that: 

•  the  sea  scatter  Doppler  spectrum  for  open-ocean  conditions  is  probably  double  peaked  for 
both  HH  and  VV  polarizations; 

•  the  ratio  of  the  Doppler  shifts  is  1.4; 

•  the  larger  Doppler  shifted  peak  (dominant  for  HH  polarization)  is  associated  with  scatterers 
near  the  crest  of  the  wave; 

•  for  X-band,  the  largest  amplitude  spectral  peak  is  the  larger  Doppler  shifted  of  the  pair  for 
HH  polarization,  and  the  smaller  Doppler  shifted  for  VV  polarization; 

•  for  lower  radar  frequencies,  the  higher  Doppler-shifted  return  may  never  dominate  the  spec¬ 
trum  for  either  HH  or  VV  polarization  because  the  scattering  mechanism  from  crests  at  these 
frequencies  is  not  as  that  strong  as  that  at  X-band; 

•  the  two  types  of  scattering  mechanisms  can  be  differentiated  by  statistical  distributions  or 
Doppler  spectra  for  the  two  polarizations 

•  with  such  sorting,  the  statistics  of  the  single-pulse  returns  of  the  high-amplitude  or  sea-spike 
fraction  of  all  returns  can  be  associated  with  the  larger  Doppler  shift  in  the  spectrum  produced 
by  coherently  integrating  the  pulse  samples  over  many  long  waves. 

Therefore,  considering  Doppler  spectra  together  with  dual  polarized  NRCS  statistics  provides  a 
preliminary  model  for  scattering  mechanisms  that  are  in  joint  agreement  and  which  together  provide 
more  information  than  each  considered  alone.  These  results  aid  in  the  interpretation  of  the  results 
presented  here  as  well.  Although  the  different  scattering  mechanisms  provide  different  statistics  for 
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the  two  polarizations  for  the  high  spatial  resolution  case  considered  here,  the  transition  to  a  Rayleigh 
distribution  with  increasing  radar  cell  size  by  the  central  limit  theorem  is  not  expected  to  be  easily 

modeled. 

SUMMARY 

We  have  reported  on  open  ocean  measurements  of  low  grazing  angle  radar  sea  scatter  made  by 
using  a  shipboard  marine  radar.  The  data  analysis  consists  of  assembling  cumulative  distributions  of 
NRCS  values  for  each  of  several  range  bins  from  an  azimuthal  sector  60°  wide.  This  sector  is  cen¬ 
tered  about  the  direction  that  exhibits  the  maximum  scatter  amplitude.  The  cumulative  distributions 
are  plotted  on  Rayleigh  probability  paper  on  which  Weibull  distributions  plot  as  straight  lines,  the 
Rayleigh  distribution  being  one  of  the  family  of  Weibull  distributions. 

The  distributions  exhibit  linear  features  in  the  high  and  middle  regions  of  the  NRCS-percentage 
plane,  while  the  lowest  amplitude  region  round  off  to  a  minimum  NRCS  value.  This  latter  rounded 
contribution  is  attributed  to  a  partial  contribution  of  noise  values  that  are  the  result  of  samples  col¬ 
lected  from  shadowed  regions  of  the  surface  and/or  samples  from  regions  that  exhibit  NRCS  values 
that  fall  below  the  noise  level  of  the  receiver.  The  two  regions  showing  linear  behavior  are  attributed 
to  two  different  scattering  mechanisms.  From  the  behavior  of  the  data  with  increasing  range,  these 
are  shown  to  be  due  to  a  distributed  scattering  mechanism  for  the  midlevel  NRCS  region  of  scatter 
and  localized  discrete  scattering  surface  features  for  the  highest  level  region  of  NRCS  values. 
Current  models  for  scatter  suggest  that  these  two  mechanisms  are  Bragg  and  sea  spike  scattering 
mechanisms,  respectively. 

The  characterization  of  the  two  distributions  with  environmental  parameters  is  accomplished  by 
determining  parameters  of  the  Weibull  linear  fits  to  the  data  and  plotting  these  vs  wind  speed  meas¬ 
ured  aboard  the  ship.  The  parameters  used  are  the  median  value  and  slope  of  the  Bragg-related 
Weibull  distribution,  and  “effective  mean”  and  slope  of  the  sea  spike-related  Weibull  distribution. 
The  intercept  of  the  two  Weibull  linear  fits  has  also  been  calculated,  and  the  NRCS  and  percentage 
value  coordinates  of  this  point  are  stored.  The  Bragg  median  value  typically  corresponds  to  the 
median  value  of  the  entire  distribution  as  well. 

The  results  of  the  correlation  of  radar  parameters  against  environmental  parameters  is  summar¬ 
ized  as  follows: 

•  The  median  NRCS  of  the  distributed  scatter  contribution  to  the  cumulative  distribution  can 
be  fit  as  exponentially  varying  with  wind  speed,  with  a  power  ranging  from  2  to  5.75  as 
the  depression  angle  varies  from  1°  to  7.7°. 

•  The  median  NRCS  varies  very  little  with  increasing  radar  depression  angle,  in  agreement 
with  a  composite  scatter  model  incorporating  Bragg  scatter  in  this  region. 

•  The  percentage  occurrence  of  sea  spikes  can  be  roughly  fitted  by  an  exponential  wind 
speed  asymptote  satisfying 

S5%  =  70%  *(t/ / U  10m/s)3  75, 


the  power  law  dependence  similar  to  that  observed  for  whitecap  coverage  of  the  ocean  sur¬ 
face.  For  the  higher  winds  in  these  observations,  the  data  fall  below  the  asymptote,  caused 
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either  by  the  seas  not  being  fully  developed  for  lack  of  fetch  or  time  or  because  the  linear 
curve  is  an  incorrect  form  in  this  region. 

•  For  all  depression  angles  >2°,  the  Weibull  shape  parameter  B  for  the  distributed  scatter  com¬ 
ponent  at  10  m/s  wind  speed  is  ~3,  increasing  ’vith  wind  speed. 

•  For  all  depression  angles  >2.5°,  the  Weibull  shape  parameter  B  for  the  sea-spike  scatter 
component  at  10-tn/s  wind  speed  is  ~2,  also  increasing  with  wind  speed. 

•  For  depression  angles  <2°,  the  Weibull  shape  parameters  appear  to  approach  unity  for  the 
distributed  scatter  component  (Rayleigh  distribution  of  radar  echo  field  strength)  and  0.75  for 
the  sea-spike  component. 

•  The  weak  variation  of  distributed  scatter  NRCS  with  grazing  angle  is  in  agreement  with  a 
Bragg  scatter  mechanism  if  tilt  modulation  is  assumed  and  by  using  local  long  wave  slope  of 
the  order  of  16°  to  20°  on  the  leeward  slope  of  the  wave.  However,  the  levels  are  much  too 
high  to  be  explained  by  horizontally  polarized  Bragg  scatter  composite  model. 

•  The  variation  of  mean  sea-spike  absolute  RCS  has  the  wrong  sense  with  depression  angle  to 
be  described  by  Bragg  scatter,  since  it  is  negative  for  increasing  depression  angle,  varying  as 
0  to  2.33  for  depression  angles  >2°.  At  this  point,  shadowing  effects  modify  the  distribution 
of  NRCS  with  decreasing  depression  angle  and  affect  the  calculation  of  the  mean  RCS.  This 
variation  with  depression  angle  is  in  the  opposite  sense  of  that  predicted  by  either  the  wedge 
model  or  Wetzel’s  spilling  plume  model. 

•  The  variation  of  distributed-scatter  NRCS  with  wind-friction  velocity,  which  attempts  to 
account  for  atmospheric  stability  effects  by  an  air-sea  temperature-dependent  drag  coefficient, 
shows  too  much  scatter  to  differentiate  between  three  current  temperature-dependent  models. 

Some  statements  can  be  made  about  the  errois  involved  ilic  Jam  analysis.  The  accuracy  with 
which  the  wind  speed  can  be  measured  is  a  topic  of  concern,  and  several  errors  can  occur.  These  are 
primarily  caused  by  the  effects  of  flow  distortion  of  the  wind  around  the  ship’s  superstructure, 
enhanced  for  cases  when  the  ship  is  not  traveling  into  the  wind.  Recent  model  studies  by  Blanc  [29] 
have  shown  that  flow  distortion  can  produce  errors  as  large  as  40%,  depending  on  where  the 
anemometer  is  located  relative  to  the  superstructure.  Since  the  relative  wind  direction  was  not  avail¬ 
able,  this  effect  could  not  be  minimized.  Future  experiments  are  planned  in  which  the  relative  wind 
vector  will  be  measured,  i.e.,  the  ship’s  heading  as  well  as  wind  direction. 

Current  models  for  high  grazing  angle  scatterometry  suggest  that  the  radar  sea  echo  may  be 
more  closely  related  to  the  wind  stress  or  wind-friction  velocity.  This  is  most  accurately  determined 
by  measuring  the  fluctuations  in  the  wind  with  in  situ  sensors  and  is  independent  of  the  absolute 
values  of  the  wind  speed.  We  have  used  an  empirical  bulk  method  of  determining  wind  friction  velo¬ 
city  in  this  work,  based  on  measurement  of  the  wind  speed.  While  wind  fluctuation  measurements 
are  also  affected  by  ship  flow  distortion,  the  errors  can  be  minimized  by  sorting  according  to  wind 
direction  relative  to  the  ship  heading,  and  these  experiments  are  also  planned.  They  will  use  dissipa¬ 
tion  techniques  to  measure  wind  stress. 

The  errors  in  the  absolute  measure  of  radar  cross  section  are  somewhat  uncertain,  since  the  cali¬ 
bration  sphere  shows  some  fluctuations,  presumably  caused  by  multipath  from  the  sea  surface.  How¬ 
ever,  preliminary  processing  of  a  second  experiment  using  S-  and  X-band  radars  on  the  same  ship 
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show  differences  of  <2  dB  for  mild  winds  when  the  Bragg  scatter  model  is  expected  to  hold.  In 
addition,  these  results  fit  the  straight  lines  determined  from  the  data  reported  here  quite  well.  Thus,  it 
is  estimated  that  the  radar  cross  sections  are  accurate  to  within  less  than  a  few  dB.  In  future  experi¬ 
ments,  upwind-downwind  and  upwind-crosswind  results  will  be  available  as  well. 
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Appendix  A 

LOW  GRAZING  ANGLE  SEA  SCATTER  BACKGROUND 


COMPLEXITIES  ENCOUNTERED  FOR  SHALLOW  GRAZING  ANGLES 

In  this  section  we  discuss  the  difficulties  encountered  for  low  grazing  angle  sea  scatter,  as  com¬ 
pared  to  that  for  higher  illumination  angles  appropriate  to  aircraft  or  satellite  radars.  Assume  that  a 
distributed  scattering  mechanism  is  operative  in  the  low  grazing  angle  case,  such  as  Bragg  scatter, 
since  it  appears  to  give  a  qualitatively  correct  description  for  higher  angles.  The  behavior  predicted 
by  this  mechanism  for  low  grazing  conditions  is  examined.  In  addition,  two  other  items  peculiar  to 
low  grazing  angles  arise:  shadowing  of  the  surface  by  crests  of  the  steepest  waves  and  the  occurrence 
of  radar  sea  spikes.  These  latter  are  deterministically  associated  with  breaking  or  sharply  crested 
features  on  the  surface  and  are  expected  to  be  the  result  of  a  localized  scattering  mechanism.  Each  of 
these  is  now  examined  in  some  detail. 

Bragg  Scatter 

At  very  low  grazing  angles,  say  <10°,  the  variation  of  the  normalized  radar  cross  section 
(NRCS)  with  depression  angle  predicted  by  the  Bragg  scatter  model  is  more  severe  than  for  moderate 
angles,  roughly  between  30°  and  60°.  This  behavior  is  shown  in  Fig.  A1  (from  Valenzuela  [A  1  ])  for 
vertical  polarization.  For  very  high  depression  angles,  between  60°  and  90°,  specular  reflections 
from  small-scale  facets  on  the  long-wave  features  of  the  ocean  surface  dominate  the  scatter.  For 
smaller  depression  angles,  in  the  60°  region  and  lower,  Bragg  scatter  is  predicted  as  dominating  the 
radar  return.  Returns  expected  from  sea  spikes  are  not  shown  here.  The  Bragg  scattering  mechanism 
is  caused  by  resonant-type  reflections  from  ocean  wave  components  with  wavelength  satisfying  the 
Bragg  relation  between  the  radar  wavenumber  k  and  the  ocean  wavenumber  K : 

K  =  2k  sin  0  (Al) 

where  0  is  the  incidence  angle  as  measured  from  nadir.  It  was  first  suggested  as  a  mechanism  for 
radar  scatter  by  Crombie,  who  was  studying  sea  scatter  at  HF  frequencies  [A2].  For  our  purposes, 
the  Bragg  model  is  sufficient  to  illustrate  the  effects  of  the  rapid  decrease  of  radar  echo  with  decreas¬ 
ing  depression  angles. 

The  NRCS  of  the  sea  surface  for  Bragg  scatter  [A3]  is  given  by 

<j°  =  47r/t4  cos4  0  I  Gjj  |2  W(2k  sin  0,0),  (A2) 


where 
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Fig.  A1  —  The  Bragg  scatter  model  is  shown  as  the  region  from  80°  downward  to  graz¬ 
ing,  in  this  plot  showing  the  variation  of  NRCS  with  local  incidence  angle  at  the  surface 


and  where  W(kx,ky)  is  the  two-dimensional  Cartesian  ocean  wave  number  spectral  density.  The 
strong  variation  of  Eq.  (B2)  at  small  depression  angles  is  also  shown  in  Fig.  Al.  However,  rather 
than  the  depression  angle  as  measured  at  the  radar,  the  local  incidence  angle  at  the  point  of  scattering 
is  important.  For  the  open  ocean  and  at  microwave  frequencies,  shallow  local  grazing  angles  are 
encountered  infrequently  because  of  the  tilting  of  small-scale  structure  by  long  waves.  Such  tilting 
will  modify  the  NRCS  observed  at  shallow  grazing,  and  this  predicted  behavior  is  discussed  in  the 
composite  model  section. 

Shadowing  of  Radar  Energy  by  Steep  Wave  Crests 

Regarding  shadowing,  radar  energy  propagating  just  over  the  tops  of  the  highest  waves  will  not 
illuminate  the  ocean  surface  just  behind  the  wave  crests  for  very  small  depression  angles.  This  effect 
is  shown  in  Fig.  A2(a-c)  for  three  different  depression  angles  and  illumination  in  the  upwind  direc¬ 
tion.  In  the  first  case,  appropriate  to  illumination  from  an  aircraft  or  satellite  radar,  the  entire  surface 
is  illuminated  and  no  shadowing  is  encountered.  In  the  second  example,  some  shadowing  of  the 
troughs  of  the  waves  occurs,  but  virtually  all  wave  crests  on  the  forward  face  of  the  long  wave  are 
illuminated.  In  the  last  example,  only  the  highest  wave  crests  are  illuminated,  these  being  associated 
with  the  dominant  waves  of  the  ocean  wave  spectrum.  The  transition  from  case  two  to  case  three  will 
be  important  for  later  considerations.  The  shaded  region  in  each  of  the  figures  is  defined  as  the  sha¬ 
dowed  region,  and  radar  samples  at  the  time  delays  associated  with  these  regions  would  result  in  sam¬ 
ples  of  receiver  noise,  since  no  echo  from  the  surface  occurs.  For  a  long-pulse  radar  operating  at  low 
grazing  angles,  such  shadowing  prevents  one  from  relating  the  NRCS  to  the  roughness  of  the  ocean 
because  the  percentage  of  the  ocean  surface  illuminated  is  uncertain. 
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Fig.  A2  —  (a)  The  geometry  for  moderate  grazing  angles  and  a  high-resolution  radar. 
Small-scale,  wind-driven,  and  surface  turbulence  fine  structure  varies  across  the  long 
wave,  producing  a  distributed  scatter  source.  Surface  features  due  to  incipient  breaking 
and  wave  crests  provide  localized  scatterers,  which  are  not  important  at  moderate  to 
high  angles,  (b)  The  same  surface  is  shown  illuminated  at  a  moderately  low  grazing 
angle.  The  surface  below  the  hatched  region  is  shadowed  by  a  wave  crest  and  produces 
samples  of  the  receiver  noise.  At  this  angle,  the  crests  and  leeward-face  breaking 
features  are  responsible  for  sea  spike  returns,  none  of  which  are  shadowed  for  this  graz¬ 
ing  angle,  (c)  The  same  surface  is  shown  illuminated  for  a  very  low  grazing  angle. 
Now  the  smaller  wave  crests  on  the  leeward  face  are  shadowed,  in  addition  to  a  large 
portion  of  the  distributed  scatter  region,  and  the  surface  behaves  in  a  manner  similar  to 
Wetzel’s  scattering  model  [20]. 


Sea  Spikes 

Radar  sea  spikes  have  been  defined  by  their  behavior  in  a  sea  echo  time  series  as  the  highest 
amplitude  returns  whose  temporal  behavior  appears  to  have  a  very  small  variation  over  an  extended 
period  of  time.  However,  as  they  are  the  largest  echo  of  the  time  series,  they  have  also  been  shown 
to  have  a  different  statistical  characterization  when  displayed  on  a  cumulative  distribution  of  RCS  by 
Hansen  and  Cavaleri  [A4]  and  by  Olin  [A5].  They  were  shown  to  be  deterministically  associated 
with  crest  and  breaking  features  on  the  surface  [A4],  and  the  scattering  mechanism  describing  them 
would  have  to  be  a  localized  one.  Two  models  exist  in  the  literature,  the  spilling  breaker  model  by 
Wetzel  [20],  and  any  one  of  many  wedge  models,  the  most  recent  that  of  Lyzenga  et  al.  [  A6] .  In  this 
report  we  have  suggested  a  modified  version  of  Wetzel's  to  explain  an  anomoly  in  much  of  the  data  at 
2°  grazing  angle. 

THE  COMPOSITE  SCATTER  MODEL  FOR  SHALLOW  GRAZING  ANGLES 

For  cases  when  waves  much  longer  than  the  Bragg  scatter  wavelength  are  present  along  with  the 
small-scale  structure,  such  as  swell  in  the  open  ocean  or  a  paddle-generated  wave  in  a  wave  tank 
experiment,  the  composite  or  two-scale  model  is  necessary  to  describe  the  radar  echo  [ A7] .  The 
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composite  model  has  been  assumed  to  hold  for  complex  ocean  waves  as  well,  and  is  the  basis  of 
attempts  to  invert  synthetic  aperture  radar  (SAR)  sea  surface  imagery  to  obtain  directional  wave  spec¬ 
tral  estimates.  As  pointed  out  by  Wetzel  [A8] ,  however,  it  is  not  obvious  that  the  more  complicated 
small-scale  structure  developed  in  the  open  ocean  can  be  described  by  Bragg  scatter.  Nonlinear 
features  generated  on  the  surface  do  not  satisfy  the  mathematical  constraints  of  the  Rayleigh 
hypothesis,  ka  <  1,  where  k  is  the  radar  wave  number  and  a  is  the  sinusoidal  component  amplitude. 
Nonetheless,  we  invoke  the  Bragg-scatter  composite  model  later  for  lack  of  another,  to  describe  the 
scatter  from  the  distributed  roughness  that  lies  along  the  greatest  fractional  area  of  the  wave,  away 
from  the  sharply  defined  wave  crests.  The  crest  scatter  may  be  less  important  for  moderate  depres¬ 
sion  angles  than  for  near  grazing  angles  and  probably  contributes  little  from  high-resolution  radars 
because  of  the  small  surface  area  of  such  scatterers.  For  a  high-resolution  radar,  however,  breaking 
wave  crests  may  cover  a  large  fraction  of  a  radar  cell  and  contribute  a  significant  percentage  to  the 
statistics  of  echo  amplitudes.  For  example,  <r°  value  of  —30  dB  in  the  10°  depression  angle  region  of 
Fig.  A1  results  in  a  0. 1-m2  cross  section  for  a  typical  100-m2  scattering  area  for  a  marine  radar.  A 
single  breaking  crest  can  produce  a  10-m2  echo  quite  easily.  Even  for  a  much  larger  scattering  patch, 
it  is  easy  to  see  how  crest  features  can  contribute  significantly  to  the  total  cross  section. 

The  composite  radar  scatter  model  includes  Bragg  scatter  as  an  element  but  allows  this  scatter 
from  the  small-scale  structure  to  vary  along  the  long-wave  profile.  This  is  caused  by  two  mechan¬ 
isms:  geometric  or  tilt  modulation,  which  simply  describes  the  local  small-scale  structure  being  tilted 
by  the  long  wave  profile;  and  hydrodynamic  modulation,  which  includes  all  fluid  dynamic  effects  that 
cause  the  spectral  density  of  the  small-scale  structure  to  vary  along  the  long-wave  profile.  To  this 
traditional  composite  scatter  model,  to  accommodate  our  interests  in  low-grazing  angle  scatter,  we  add 
the  broad  range  of  effects  that  we  label  chaotic.  These  effects  include  wave  breaking,  which  is  par¬ 
ticularly  important  near  and  just  forward  of  the  crests.  This  mechanism  has  not  been  treated  by  pre¬ 
vious  composite  models  but  is  necessary  to  describe  the  occurrence  of  sea  spikes  found  in  low  grazing 
angle  radar  sea  scatter. 

Geometric  or  Tilt  Modulation 

Whereas  the  expression  in  Eq.  (Al)  holds  for  each  wave  component  moving  on  a  locally  flat 
surface,  the  sea  surface  is  composed  of  many  wavelengths  and  heights.  As  a  result,  the  local 
incidence  angle  at  the  surface  is  not  the  same  as  the  depression  angle  at  the  radar  but  is  modified  by 
the  local  slope  of  the  long  wave  at  the  point  of  radar  illumination,  as  is  seen  in  Fig.  A3.  Thus,  the 
expression  for  the  Bragg  scatter  amplitude  given  in  Eq.  (Al)  requires  this  local  incidence  angle, 
which  is  the  sum  of  the  grazing  angle  and  the  slope  angle  of  the  long  wave  component,  with  slope 
measured  positive  as  defined  in  Fig.  A3. 


TRUE 


Fig.  A3  —  Local  incidence  angle  of  the  radar  signal  is  determined  by  the 
depression  angle  at  the  radar  and  the  local  slope  of  the  long  wave  at  the  point  of 
incidence  and  causes  tilt  modulation  of  the  expression  of  Eq.  (3) 
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For  a  collection  of  samples  at  a  fixed  delay  and  azimuth  for  a  high  spatial  resolution  radar,  a 
large  range  of  echo  amplitudes  results  from  tilt  modulation.  The  NRCS  samples  forms  a  distribution 
of  values  given  by  Eq.  (2)  evaluated  at  angles  between  the  sum  of  the  depression  angle  at  the  radar 
plus  and  minus  the  maximum  slope  of  the  long  wave  illuminated.  As  an  example,  for  a  7-s  ocean 
wave  of  75-m  wavelength,  the  maximum  long  wave  slope  is  arctan  [2  ir(h  / 2)  / /]  =  5°,  in  effect  the 
half-angle  of  the  tilt  modulation.  Thus,  for  a  grazing  angle  in  the  plateau  region  at  45°,  local  depres¬ 
sion  angles  will  vary  between  40°  and  50°.  Using  Fig.  A3,  the  NRCS  is  modulated  by  this  geometric 
effect  between  values  of  -18  and  -15  dBm2.  For  a  grazing  angle  of  8°,  the  distribution  of  local 
grazing  angles  varies  between  3°  and  13°,  resulting  in  a  modulation  of  the  NRCS  of  Fig.  A1  of 
between  —45  and  —25  dBm2,  a  significant  increase  in  dynamic  range.  For  very  low  grazing  angles, 
the  local  incidence  angle  is  varied  from  0  (negative  local  incidence  angles  are  cut  off  due  to  the  onset 
of  shadowing  at  0°)  to  just  above  5°.  Because  the  NRCS  plunges  so  strongly  as  0°  is  approached, 
the  receiver  noise  level  is  reached  and  limits  the  smallest  NRCS  observed.  Nontheless,  the  variation 
in  NRCS  is  largest  at  the  lowest  radar  depression  angles. 

Since  the  variation  of  the  local  incidence  angle  along  the  long  wavelength  components  of  the 
ocean  surface  does  not  produce  a  severe  variation  in  amplitude  dynamic  range  for  moderate  incidence 
angles,  a  SAR  image  of  the  ocean  surface  can  be  recorded  by  using  just  a  few  bits  of  amplitude  infor¬ 
mation.  The  marine  radar,  on  the  other  hand,  typically  has  ~60  dB  of  dynamic  range  to  handle  the 
wide  range  of  echo  power  expected  at  low  grazing  angles. 

Hydrodynamic  Modulation 

In  addition  to  the  geometric  modulation,  a  hydrodynamic  modulation  also  occurs.  It  is  well 
known  that  the  small-scale  surface  spectral  density  is  modulated  by  the  orbital  wave  motions  of  the 
long  waves  [A8].  Straining  of  small  waves  by  long  waves  causes  a  shortening  of  the  wavelength  of 
the  front  face  and  a  lengthening  on  the  back  face.  A  second  type  of  modulation  is  that  caused  by 
sheltering  of  the  wind  by  the  wave  profile,  causing  a  variation  in  the  local  wind  speed  across  the  sur¬ 
face,  with  an  expected  smaller  Bragg  scatter  density  excited  in  the  sheltered  regions.  The  effect  of 
each  of  these  is  to  vary  the  amplitude  of  the  Bragg  scatterer  spectral  density  of  Eq.  (A2)  as  a  function 
of  position  on  the  long  wave.  These  effects  are  typically  dealt  with  empirically  by  definition  of  the 
modulation  transfer  function,  introduced  by  Wright  [A7].  The  modulation  transfer  function  is  defined 
as 


a(x)  =  a°[l  +  m(f)  i)(x ) ] ,  (A3) 

where  <j(x)  is  the  spatially  varying  NRCS  of  the  sea-surface  height,  jj(x).  It  consists  of  a  dc  term,  a0, 
and  a  spatially  varying  contribution  that  is  proportional  to  the  local  sea-surface  long-wave  profile. 
This  empirical  expression  assumes  that  the  mapping  from  radar  image  to  wave  amplitude  is  linear  for 
a  given  wave  frequency,  which  has  been  shown  to  be  true  in  wavetank  experiments  for  moderate 
grazing  angles.  Such  an  empirical  formalism  does  not  hold  for  low  grazing  angles  and  open-ocean 
waves  for  reasons  discussed  in  the  text. 

Chaotic  Events,  Including  Wave  Breaking 

Recent  results  have  shown  that  deterministic  features  on  the  surface  may  also  contribute  signifi¬ 
cantly  to  the  NRCS  of  the  sea  surface  for  the  higher  grazing  angles  in  which  specular  scatter  may 
occur  [A9],  particularly  for  the  very-high-resolution  conditions  encountered  in  wavetank  studies  in 
which  a  sharp  surface  crest  feature  can  fill  the  illumination  cell.  For  low  grazing  angles  encountered 
by  shipboard  radars,  scattering  features  are  not  expected  to  present  specular  scatter,  but  a  type  of  dou¬ 
ble  scatter  can  occur  from  plume-type  breakers  rolling  down  the  front  face  of  a  propagating  ocean 
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wave  [A8].  The  effect  is  not  expected  to  be  as  dominant  as  in  a  wavetank,  however,  because  the 

feature  does  not  fill  as  large  a  fraction  of  the  scattering  cell  in  range  or  azimuth. 
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Appendix  B 

AIR-SEA  INTERACTION  PARAMETERIZATION 


WIND  STRESS  -  WIND  SPEED  -  FRICTION  VELOCITY  RELATIONSHIPS 

To  interpret  the  apparent  trends  in  the  plots  ot  radar  parameters  vs  wind  speed,  it  is  useful  to 
introduce  the  wind  stress  r  [B I] : 


T  =  p  <uw> ,  (Bl) 

where  p  is  the  ratio  of  the  air  density  to  water  density  and  <uw  >  represents  the  correlation  of  the 
two  horizontal  components  of  the  vector  wind  with  the  coordinate  system  taken  as  the  x-axis  along 
the  wind  direction.  The  wind  friction  velocity  u*  is  then  defined  by 

u*  =  <uw>  (B2) 

for  purposes  of  bulk  or  large-area  average  parameterization  of  the  air-sea  interaction.  The  wind  fric¬ 
tion  velocity  is  the  large-area  averaged  value  of  the  wind  speed  at  the  skin  of  the  hypothetical  average 
rough  surface,  defined  so  that  the  stress  is  uniformly  distributed  over  that  surface.  The  physical  real¬ 
ity  of  such  a  quantity  on  a  small  scale  is,  of  course,  not  meaningful  in  the  case  of  a  two-scale  sea  sur¬ 
face  from  which  radar  returns  occur.  Nevertheless,  a  formal  definition  can  be  made,  considering  the 
stress  and  friction  velocity  averaged  over  an  area  containing  many  large-scale  wavelengths. 

The  friction  velocity  can  be  more  effectively  related  to  the  normalized  radar  cross  section 
(NRCS)  cumulative  distributions  than  to  the  wind  speed  for  two  reasons.  First,  the  wind  stress  defin¬ 
ition  includes  the  effects  of  atmospheric  stability  or  air-sea  temperature  differences.  Second,  the  con¬ 
tributions  to  stress  delivered  from  the  wind  to  the  surface  can  be  broken  up  into  several  contributions, 
each  of  which  may  be  able  to  be  related  to  different  contributions  to  the  NRCS  cumulative  distribu¬ 
tions.  These  wind  stress  contributions  serve  as  some  of  the  source  and  sink  terms  in  the  equation  for 
wind  wave  generation  and  energy  transfer  between  wave  components  of  the  wave  spectrum  [B2], 
Intuitively,  it  would  not  be  surprising  for  a  relationship  to  exist  between  sea  surface  features/radar 
scattering  mechanisms  and  individual  surface  stress  contributions. 

AIR-SEA  TEMPERATURE  DIFFERENCE  EFFECTS 

Air-sea  temperature  differences  have  been  reported  as  an  important  parameter  in  the  measure  of 
mean  NRCS  for  single  wind  speed  [B3] .  The  effect  of  air-sea  temperature  difference  in  the  wind 
stress  definition  can  be  taken  into  account  by  a  drag  coefficient  CD  in  a  bulk  parameterization  of  the 
wind  friction  velocity.  The  bulk  method  of  measuring  the  wind  stress  assumes  that  the  wind  friction 
velocity  defined  at  the  surface  can  be  related  to  the  wind  speed  at  an  arbitrary  height  H  through  the 
drag  coefficient 
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where  U{H)  is  the  magnitude  of  the  wind  speed  measured  at  an  altitude  H ,  and  CD  is  determined 
empirically.  The  drag  coefficient  is  generally  defined  only  for  neutrally  stable  atmospheric  condi¬ 
tions;  for  moderate  winds  it  best  fits  an  equation  of  the  form 


CD  =  a  +  bU.  (B4) 

Definition  of  the  drag  coefficient  for  unstable  atmospheric  conditions  has  only  recently  been 
done,  and  at  least  four  different  formulations  are  available  for  comparison  [B3-B6],  If  such  a  defini¬ 
tion  can  be  made  successfully  and  if  the  spread  in  the  radar  parameters  previously  discussed  is  caused 
by  the  air-sea  temperature  difference  dependence  of  the  drag  coefficient  alone,  the  spread  of  the  radar 
data  is  reduced  significantly  when  plotted  against  wind  friction  velocity  compared  to  the  wind  speed. 
The  definition  made  by  Keller  et  al.  [B3]  is  not  successful  in  reducing  the  spread  of  the  mean  RCS 
for  X-band  radar  scatter  at  45°  incidence  for  measurements  made  in  the  Gulf  of  Mexico. 

WIND  STRESS  AS  THE  SUM  OF  SEVERAL  SEPARABLE  CONTRIBUTIONS 

Donelan  [B7]  has  suggested  that  while  the  wind  stress  is  a  highly  averaged  quantity  as  formally 
defined,  one  can  think  of  it  as  consisting  of  a  sum  of  terms,  each  associated  with  different  physical 
air-sea  interaction  processes.  These  are  first  broken  into  viscous  drag  ts  and  form  drag  contribu¬ 
tions.  The  form  drag  is  broken  into  terms  related  to  production  of  surface  currents  rc  and  production 
of  momentum  in  the  form  of  surface  waves.  The  momentum  term  further  consists  of  a  source  term 
tm  and  a  sink  term  tb  relating  to  momentum  lost  through  wave  breaking.  The  final  expression  is 
simply  the  sum  of 


tto,  =  rs  +  tc  +  tm  +  tb.  (B5) 

In  particular,  he  has  shown  by  simultaneous  measurements  of  wave  breaking  and  of  the  fluctuations  of 
the  water  particle  motions  below  the  surface  <  u'w'  >  that  breaking  waves  deliver  a  significant  frac¬ 
tion  of  momentum  flux  through  the  surface  into  the  water,  resulting  in  underwater  turbulence.  Phil¬ 
lips  [B2]  has  further  suggested  that  the  contribution  to  the  wind  stress  due  to  wave  breaking  can  be  a 
significant  one  under  fully  developed  seas  based  upon  theoretical  considerations.  We  continue  with 
this  wavebreaking  contribution  to  surface  stress,  following  up  on  an  observation  by  Jin  Wu  [B8]. 

Wu  reanalyzed  existing  data  in  the  literature  on  the  wind  speed  dependence  of  percentage-area 
of  whitecap  coverage  on  the  ocean  surface  and  inferentially  argued  the  reason  for  a  dependence  of  the 
drag  coefficient  as  wind  speed  to  the  one-half  power  to  fit  a  3.75  slope  of  log  of  whitecap  coverage  to 
wind  speed.  The  argument  was  that  the  whitecap  coverage  should  be  proportional  to  the  energy  flux 
into  the  surface  under  equilibrium  conditions.  This  energy  flux  is  the  product  of  the  wind  stress  times 
the  wind  drift  currents,  the  wind  drift  currents  being  also  proportional  to  the  wind  friction  velocity. 
Wu  thus  proposed  the  following  dependence  for  the  percentage  area  coverage  by  whitecaps  W: 

W  ~  tu.  =  pn,3  ~  (CDU2)y2.  (B6) 

If  one  further  assumes  a  wind  speed  dependence  of  the  drag  coefficient  of  U  to  the  first  power,  as 
found  by  Large  and  Pond  [B9J,  then  one  finds  the  energy  flux  and  W  proportional  to  U45.  To 
satisfy  the  observed  whitecap  coverage  dependence  of  f/3  75,  Wu  had  to  postulate  a  wind  friction  velo¬ 
city  dependence  on  wind  speed  of  U0  5.  Such  a  dependence  agrees  neither  with  the  measurements  of 
Large  and  Pond  [B9]  and  others,  nor  with  previously  published  work  by  Wu  [BIO], 


NRL  REPORT  9099 


REFERENCES 

Bl.  H.U.  Roll,  Physics  of  the  Marine  Atmosphere  (Academic  Press,  New  York,  1965). 

B2.  O.M.  Phillips,  “Spectral  and  Statistical  Properties  of  the  Equilibrium  Range  in  Wind-Generated 
Waves,  I.  Fluid  Mech.  156,  505-531  (1985). 

B3.  Keller,  W.C.,  W.J.  Plant,  and  D.E.  Weissman,  “The  Dependence  of  X-band  Microwave  Sea 
Return  on  Atmospheric  Stability  and  Sea  State,”  J.  Geophys.  Res.,  90  1019-1029  (1985). 

B4.  L.B.  Wetzel,  “On  Microwave  Scattering  by  Breaking  Waves,”  in  Wave  Dynamics  and  Radio 
Probing  of  the  Ocean  Surface,  O.M.  Phillips  and  K.  Hasselmann,  ed.  (Plenum  Press,  1986  pp. 
273-284. 

B5.  D.R.  Lyzenga,  A.L.  Maffett,  and  R.A.  Shuychman,  “The  Contribution  of  Wedge  Scattering  to 
the  Radar  Cross  Section  of  the  Ocean  Surface,  IEEE  Trans.  Geosci.  Remote  Sensing,  GE-21, 
502-505  (1983). 

B6.  Jin  Wu,  “Stability  Parameters  and  Wind  Stress  Coefficients  under  Various  Atmospheric  Condi¬ 
tions,”  J.  Atmos.  Oceanic  Tech.,  (in  press),  1986. 

B7.  M.A.  Donlen,  Whitecaps  and  Momentum  Transfer,  Turbulent  Fluxes  through  the  Sea  Surface, 
Wave  Dynamics,  and  Prediction,  (Plenum  Press,  New  York,  1978)  pp.  273-287. 

B8.  Jin  Wu,  “Oceanic  Whitecaps  and  Sea  State,  J.  Phys.  Oceanog.,  9,  pp.  1064-1068  (1979). 

B9.  W.G.  Large  and  S.  Pond,  “Open  Ocean  Momentum  Flux  Measurements  in  Moderate  to  Strong 
Winds,  J.  Phys.  Oceanog.,  11,  324-336  (1981). 


B10.  Jin  Wu,  “Wind  Stress  Coefficients  over  Sea  Surface  Near  Neutral  Conditions— A  Revisit,  J. 
Phys.  Oceanog,  10,  727-740  (1980). 


